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basis of a discrete model ‘that has gaeviowsly been shown to be the finite dif- 
ference analog of the continuum model of a plate. The equations enn ae 

_ Clude the effects of flexural and torsional stiffnesses of the beams. . 


wey 


Many problems that are of interest to engineers involve boundaries or oll 
_ that are not readily suitable for analytic solution. Such situations may : 
be resolved by numerical a that provide approximate solutions to 


_ The method of finite differences is one of the methods, — This method has" 
been used increasingly in 1 many applied fields since the advent of modern n dig- 7 
‘ital computers. There are two basic approaches in seeking approximate sole- 
= to physical problems of continuum by the method of finite ral — 
and subsequently approximate the differential equation by a system of finite . 
difference equations, The other approach is to introduce a discrete model to 
approximate | the ‘continuum, with the restriction that in the limit the discrete 


Note. —Discussion open until May 1, 1962. To extend the closing date one month, a 
s itten phase no must be filed with the Executive Secretary, ASCE. This paper is part 
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model approaches the continuum (at least in an intuitive esense). 


pr 

beam systems. Using this approach, difficulties due to discontinuities and spe- 
boundary conditions that would otherwise be encountered, are circum- 
ve 


pplication of finite differences to 
What is claimed to 4 
‘son one is the use of the approach described herein (including the use of 
fictitious deflections) that makes possible the derivation ef Ge equations for 
plate- -beam systems: (based on central differences) wherein the beams have 
Comins torsional and flexural stiffnesses. In addition, most of the equations, _ 
those pertaining to corner boundaries with beams, are presented 
= object of the paper is twofold: To demonstrate that a physical model 
can be used to great advantage in supplementing a purely mathematical ap- 7 
ti proach, thus is making possible solutions to problems involving complex bound- - 
ary conditions and constraints; and to describe | the derivation of and presen = 
- equations that can be used by engineers in the analysis of plate-beam a 
involving beams with torsional and flexural stiffnesses. 
—_ analysis of plate-beam systems is based on the ordinary theories of 


flexure of beams and plates. In addition, the neutral axes of the beams are 
assumed to coincide with the middle plane of the plate. Torsion in the beam is z 
> ‘4 assumed to be uniform between node points, and effects of warping are oe 


lected. A beam is represented by a line. 
_ Notation.—The terms adopted for use in this paper are defined where ile mf 
first ‘appear and are listed alphabetically, for convenience of reference, in the — 


Although the analysis of plate- beam systems requires finite difference re- 7 
current equations for plates, 

‘They can be found in a number of in the Recurrent e equa- b 
tions that are peculiar only to combined plate-beam systems are derived. This 
is done in some detail for one case whereas others are presented without ram- ; 

ifications except notes concerning their applications. 

_ The equations are derived for equal grid distances and are, therefore, re- 

_ stricted to square nets with straight boundaries, However, equations for rec- _ 

= nets can be ‘Similarly derived without special difficulties. For other — 


38 “A Numerical Procedure for the Analysis of Continuous Plates,” by A. . = -s _ Ang 7 
and N. M. Newmark, Proceedings, 2nd ASCE Conf. on Electronic > Computation, Septem- 

: “Analysis of Skew Slabs,” by V. P. Jensen, Bulletin No. 332, Univ. of Ulinois Engrg. 
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= _ used is an extensionof the plate model first proposed by Nathan M, Newmark,? i a 
| 
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as Central differences : are used 'throughout. The inaccuracy involved in using» 
the equations, thereiore, is a result of the error inherent in the central differ- _ 
_ Effects of Beam Stiffness. —Beams with only flexural stiffness have been 
considered before.® The effects of beams with combined flexural and torsional - 
stiffnesses will be « described, 
“Fig. 1 shows a plate : stiffened with an interior beam, tf the beam has tor-_ 
" sional stiffness there will be discontinuity at the junction between the plate and 
: that is, the curvatures of the me normal to the axis of the beam 
be different at the two faces of the beam. Referring to Fig. these cur-— 
-vatures at joint o, by the three-point formula, are as follows: a 


Wp and wg Wg are fictitious, deflections at points n and s, respectively, 7 


and A is the @iutence between node | points. The two fictitious deflections Whe 
= Wg are related by the fact that the angle of twist of the beam tiated oO can be- 


(3) 


Wn 


he curvatures given 1 moments normal to 

ee two faces of the beam that, in general, are not equal, This difference in — 

bending moments about the axis of the beam acting at joint o is 


the torsional moments in the beam at | 0, as shown in Fig. 2 2(a). | That is 


* Mos) * (Toe - Tow) = 


_ 5 “Numerical Methods in Engineering,” by M.G. Salvadori and M. ancien Prentice— 
‘Hall, Inc., Englewood Cliffs,N.J.,1952,0 
_ 6 “Moments in Simply Supported Skew I- -Beam Bridges,” ” by T. Chen, 


and N. M. Newmark, B Bulletin No. 439, Univ. of Illinois Engrg. ae Sta., 1957 
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1. STANDARD DESIGNATION OF POINTS TO 


ements ont 


siaialaaini OF JOINT ND BAR IN A PLATE- BEAM — 
SYSTEM 


— 
a 
Point O is where an equation is applied 
lg 
(a) Moments about x-axis at joint (b.) ‘Shear forces at joint O. 
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- 2Wo + Ws, 


= 


IN (“= Wse Wn - Ws Ws 


renee 


a measure of the torsional stiffness of the Substituting these expres-_ 
sions into Eq. 4 and using Eq. 3 to eliminate the fictitious deflections below 
_ the beam, g gives the equilibrium of moments about the axis of the beam in terms) 
of deflections as follows: 
2 Wn J Ws + + 


or, in pattern form 


wa =0.. 


+ + (243) ) 
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‘te 
| 
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tional requirement for the joints on the beam as given by Eq. 7. The values in- 
_ gide the dotted boxes are coefficients of the fictitious deflections at the in- : 
7 dicated points. These fictitious deflections, in general, do not coincide with 
the real deflections | at the same "points. . Eq. 4, and similar equations to be 
derived, are to be used either to eliminate the fictitious deflections or to pro- | 
vide additional equations for determining the fictitious deflections. The latter | 
procedure; is usually preferable, because the fictitious deflections are required © oe 


for finding the plate bending g moments 2 at, and normal to, the e faces of the beam, 7 > 


as well as for finding the torsional moments inthe beam. =  —— 


— .| F 

stiffness per unit width of 
= 
4 

The nreseanra anf a heam with taraion etiffnecec therefore imnnees an addi. a 
| 
— 


on The fictitious deflections can be considered as providing a means of speci- — 
_ fying the angle of twist of a beam at a given point. The use of fictitious deflec- _ “ 
tions, therefore, makes possible the use of central difference oe 
_in deriving the governing for a plate- beam beam 

Joints on or Near an Interior Beam. —Referring to Fig. 2(b), the 

of forces to middle - of the plate fora joint on an interior ‘beam 


in which Q is ee Pree concentrated load at the joint.” The other terms 


in 1 Eq. 8 shears are, , for ‘example [see Fig. -2(c)] 
+ - (Mo1 + 


ihe 


- 2w, + w 


ste 


, (10b) 


a measure of the flexural stiffness o of the beam. - Note the use of fictitious de- 
- flections in Mo1 and Mo2- The other terms in Eq. 8 can be similarly expressed | 
in terms of deflections. Substituting these | expressions into Eq. 8 and making | 


use of Eq. 3 gives the equation of equilibrium of normal forces fora joint on _ 


an interior beam, as follows: ws 
(a7. 


M. Newmark, Transactions, ASCE, Vol. 108, , 1943, p. 1161. 
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+ H) + 4H) (20 + (1 + H) 


1 does not containany fictitious deflections. However, this is not thecase , 


in the corresponding equation for the joints immediately adjacent to the beam : 


as shown in Eq. 12, in which the value inside the dotted box is the coefficient 


values of the fictitious deflections, that 
tal 


um of normal forces 


+ 4H) (10+ 6H) -(4 + 4H) 


— 
— L 
— 
Ea, 1 
— 
— 
the fictitious deflection, — 

ae | Ea. 7 applied to every joint on an interior beam provides additional equations [i 
an Exterior Beam, —Equilibri for a 
ee joints on an exterior beam and immediately adjacent to it are described, re- 
spectively, by Eqs.13and14,00 
m 
— e 


8 | 


+ 
wee 


Eq. 15 should be Ye applied to every -y joint on the exterior beam to provide the nec- - 
essary additional for determining the fictitious of 


—— on Corners with Beams. —Two types s of ¢ corners with eit at are con- 

sidered, a sharp corner and a re- entrant corner, 


For a sharp of normal forces. is given by 


2- 


this case, there are two moment to be s: satisfied, 


0 and uM 17) 


[+ 
equilibrium of moments about the axis of the beam fora joint on the exte-— 
| 
— 
| 
a 
= 


rms of deflections a are, respectively 


(2H + 


For a re- -entrant corner, ‘equilibrium of normal forces is given by 


flections at ne are generally 


-p) t-2(2H + J + 1)! w=0...(18) 
— 
— 


The two fictitious deflections at point ne in Eqs. 21 and 22 are also different. 
The one in Eq. 21 defines the angle of twist of the beam oriented in the y-— 


4 direction, whereas that in Eq. 22 defines the angle c of twist of the beam oriented | : 
thex-direction, 
Eqs. 18 and 19 and Eqs. 21 and 22 provide the necessary equations for de » 
termining the corresponding fictitious deflections. __ 
In the derivation of Eqs. 20 through 22 there are certain neliitelaiaieg. 
peculiar to a re-entrant corner that deserve further explanations, For this 
_ purpose, it it will s suffice to investigate some of the plate bending mi moments acting — 
at such Conner as follows ied 3): 


ot 8) 


‘ 
‘ 
= 
7 _while the y-curvatures are different; that is, the y-curvature of mp, is the y- — 
_ curvature of the plate at the corner while the y-curvature of mog is the same — 2 
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curvature in M3 is the same as that of the other beam, but the y- curvature is 
same as the corresponding curvature ofmoj. 


PROBLEMS 


a _ Two example problems are given to illustrate the use of the equations d de- 
—) herein, _ These problems are also intended to give some indications of 
the type of complex boundary-value problems of plate-beam systems that can 

be solved. However, extensive results and analysis of the vai variables in plate- 


beam ‘Systems | are not within the Scope of the the paper. a Bip: ie 
The application of certain ‘equations may require clarification. this 
purpose, some of the equations are explicitly written for the two problems us-_ 


ing p = in of moments at certain sections of the 


Plate bending moments about 
y-axis. shown 


JOINT AT RE- -ENTRANT CORNER_ 


?P plates are shown in Figs. 4 and ‘/. In both problems, a uniformly a 
transverse load is applied and = 0 is used, 

Ray m4 1.—One quarter of ‘the structure analyzed in the first problem is —. 
a shown in Fig. 4(a), which is a square plate with two opposite edges fixed and 

_ the remaining edges stiffened with beams. One flexural stiffness of the beams ~ 

= 20) and three different values of torsional : stiffness = 0, 10, 20)are 


Referring to Fig. 4(a), note that points 1 1 through 5 and 16 through | 20: rep- 


resent points having fictitious deflections. > 


: _ Applying Eq. 12 to point 18 and making use of Eq. 3 give _ 


2iwig - (wiz + wi7 + wig + w23) +2 + w14 + w22 + ¥ 


y 
+ * wie + w20 = 


— 
u 6 
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5. -EXAMPLE PROBLEM II 


rg 
‘4 — — 


Applying Eq. 7 to potas 23 


+ wag) + (243) Wig - wig - ( 


_ The application of the other equations, such as Eq. 11 for points 21 tga 
: 25 and Eqs. 13 and 15 to points 6 through 10, ‘needs no further explanation, — 
i: _ Example 2.—Fig. 5(a), shows one quarter of the structure analyzed in this - 


problem. _The structure is a square plate supported on four infinitely rigid 


_stiffmessesofH=J=200 
‘Taking into consideration the symmetry of the structure, Eqs, 16 and 18 

~ when applied to point 12 are straight-forward and do not require further | - com- 
mente, a both Eqs. 20 and 21 to point 65 yields, respectively, the fol- 


Wer 


2 w, +2 w: 2 w56 + +4 + (2H wag 0.75 3 


q 


a as all the exterior edges are e stiffened with beams having flexural and torsional 


3 


+ 3) 1 


- should be e emphasized | that in any problem, the number of “equations is 


equal to the number of unknown deflections including th the fictitious deflections, 


~CONCLUSIONS 
‘By introducing approximations in the physical model itself, rather than in 
the governing differential equation, finite difference recurrent equations for 
plate-beam systems can be derived without the difficulties due to discontin- _ 
uities that would otherwise be encountered, These recurrent equations when | 
used in conjunction with the equations for plates provide a means of obtaining 


_ numerical solutions to combined plate- beam systems. 


oath equations contained herein were derived on the basis of a plate model 
th 


- has ae shown to to be the finite difference equivalent of a continuum. 3 
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The following notation, adopted for use in this ‘paper, 


@6e 


_of reference and for the convenience of discussers: 


= I 

: 2: Nn? a measure of the flexural stiffness ofa a beam; 


x? ,» a measure of the torsional stiffness of abeam; 
s = bending moment in a beam; a a 

moment per given width of plate; 


— , stiffness per unit width of plate; 


= equivalent concentrated load ata Joint; . 


= = intensity of uniformly, distributed load; 
= torsional moment in a beam; 


= shear force; 


= 


= distance between joints | or node points; and - 


= Poisson’s ratio. > 
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LECTRONIC ANALOG G COMPUTERS FOR LAMINAR FLOW PR PROBLEMS 


Young! 


ivi- 


A simple techniave is presented for the use of an analog com-— 
puter to determine velocity distributions in both steady and unsteady uniform. 
laminar flow problems. Three specific examples are given, and a comparison ~ 
of the computer solutions with thet! theoretical solutions indicates the satisfac- 
tory accuracy of the technique. 


Electronic differential analyzers, or more simple ‘ “analog computers,” are 
po “now being used extensively in practically every branch of engineering, particu- 
larly for problems in which the independent variable is time. Because the 
- analog computer integrates v with ‘respect to the single variable time, it is well 
_ Suited for this type of problem. However, for many problems in areas such as . 
heat transfer, aerodynamics, and nuclear engineering, ‘solutions to 


y ‘required, Two techniques are currently available for solving partial differ- 
_ entialequations with the analog computer. In certain instances, the partial dif-— 
_ ferential equation can be reduced by separation of variables to a system of or- 
_ dinary differential equations, each involving a single independent variabl 


These equations can then be solved on the computer, However, a frequent] a 
_ Note.—Discussion open until May 1, 1962. To extend the closing date one month, a 7 
written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the — 
_ American Society of Civil Engineers, Vol. 87, No. EM 6, December, 1961. 2 
Prof. of Theoretical and Applied Mechanics, Dept. of Nuclear Towa 
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December, 1961 
simpler and more general sneniitnen te one of replacing the original partial 
derivatives by their finite -difference approximations, thereby reducing the 
original differential equation to a form that can be solved on the computer, | 
This method was first described by T. A. Rogers2 and R. M. Howe and V. S. 
Haneman3 and has for the solution of sev several different t types of of par- 


_ Although conventional electrical analogies, , such as systems utilizing con- 


the study of fluid particularly in the study of the flow of 
4 ideal fluids (see for example the work by Malavard#), no reference to studies ~ 
- in whichan electronic analog computer had been used to investigate the steady | 


or unsteady viscous flow of an incompressible fluid through or between various 
_ regular geometries such as flat plates, 


circular tubes, or rectangular tubes © 


problem is similar to the heat equation which has been studied by the use of an 


he found, Because the governing partial differential equation for this type of 


electronic analog computer, 2 laminar flow problems ; can also be treated in this 
manner, It should be noted that passive analog computers, as demonstrated by - 
J. L Shearer, 5 can also be used for this type of problem, | However, it was the ¥ 
I purpose > of the present investigation to demonstrate the applicability of the use 7 
of electronic analog computers to various types of laminar flow problems, Ih 
= this study three three specific problems —_ considered (1 (Fig. 1 ee 


(1) Starting between infinite parallel plates; 
(2) flow between infinite parallel plates under the influence of a pulsat-— 


(3) flow near an os¢ 


These three relatively simple cases were selected so that the analog solu- 
_ tions could be > compared with the theoretical solutions. _ However, as noted in a 
the last section of this paper the general n method of approach can be extended 
to through tubes of various cross sections. 
jotation.—The letter symbols adopted for use in this paper are defined 
they first appear, in or in the text, ae are arranged 


2 “Electronic Analog Computers and Partial Differential Equation Solutions,” by T, 
A ‘Rogers, Memorandum to Dept. of Engrg., Univ. of California, Los Angeles, 1952. 
"(Noted in reference 3. Memorandum no longer available) __ 
3 The Solution of Partial Differential Equations by Difference Methods Using the 
Electronic Differential Analyzer,” by R. M. Howe and V. S. Haneman, Proceedings, — 7 
“The Use of Rheoelectrical Analogies in Aerodynamics, by L. C. Malavard, AGARD 
Poa pong re Network Analog Study of Viscous F low Normal to Parallel Evenly ‘Spaced 
J. Research Journal, Vol. . 29, 1959, p. 467. 
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December, 1961 
the kinematic viscosity, u denotes a: velocity component. in in — 
‘&- direction, is time, p is fluid density, and 22 is s the pressure gradient. 
The boundary conditions for this problem are simply that fluid velocities - _ 


at the boundaries a are the same as the velocity of the boundary; that is, fora 
= Uo fo (t) as illustrated 
in Case III of Fig. 

_ To rewrite Eq. 1 in terms finite the flow 
=. field is subdivided into a number of intervals (Fig. 2). The average velocity 


gradient between, for example, nodes 4 and 5, 
j 


betwee: een 5 6 


The second derivative at at node Be can be be approximated as 


56 445] (Ay) 


Therefore, at node 5 ‘the differential equation can be ‘approximated 


um 1 ap 


ae For the particular flow problems studied eight intervals were used, For 
r Cases I and II (Fig. 1) the intervals were equally spaced, but the spacing for 7 
7 Case III was not equal because this would obviously require an infinite number» 


a and, in general, the finite difference ee 
— 
_ Eq. 6 is valid for any node point around In cer-_ a 
tain instances it is convenient to use unequal spacing, and the finite difference 
a _ equation must be modified accordingly. For the example illustrated in Fig. * alti 
2(b) the approximation for the second deriva 
(ay)? | — 
rat 


The spacing for Case III was arbitrarily selectedso that intervals” 


_ between nodes 0 and 5 were equally spaced, the interval between 5 and 6 equal 
to2Aa y, between 6 and 7 equal to 8 Ay, and between 7 and 8 equal to n Ay in 

: which 7 > >> 8. . The actual method used to simulate this condition will will be ex-— 
amined in the following section. —_ 


if Eq. 6 is rewritten as 


dp 
which ch F(t) = pee with the equa the out output voltage of 


operational of Fig. 3, the between the equations is 


in which V represents voltage, R resistance, C capacitance, @ an2 arbitrary 
= and 7 computer time. Thus it is noted that, if the correct inputs are 
: fed into the amplifier, the output can be made proportional tO Um- _ This can 
be accomplished by letting 


(ay)? F(t) 


_ then: the voltage is proportional to velocity as as indicated by Eq. 10. ~The con- 
stant of integration has been omitted from Eq. 8, because for all flow cases it = 
-. was assumed that at t=c ‘0, U=O. _ In the computer setup this corresponded to hav- 
7. no initial charge on the capacitors. The complete basic computer | circuit | 
that was used is shown in Fig. 4. The velocities correspond to the output volt- 
a ages of the gg but are the negatives of the voltages for amplifiers Al, - i 


‘dams the sign of the input voltage. If ‘desired, : sign ‘changers could be incor 
into the circuit, but this would require four additional amplifiers. ‘The 


{ 
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— 
in which ny and ny areconstants, 
Substitution of Eqs. 10, 11, and 12 into Eq. 8 indicates that 
= 
x 
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requenc’ on ge : 
computer, a Hewlett Packard low frequency 

_ Because each of the cases studied required some modification of the basic 
- computer set-up, the details of the computer circuit for each will be picusineaiinadl 


The actual computer results are given in the next section. 


(a) Case L—For this the boundary conditions are t) = 
which corresponds to Vo = — 0, and the pressure gradient is constant ‘so > 
that the input voltage, Vp, - is a constant. Because the transient flow was of 
interest for this case, the procedure followed was to apply the constant voltage - 
; Vp suddenly, and to record the output voltages of all the amplifiers continuously 
+ until a steady state condition was attained. aoe 
of (b) Case I. .—As for Case I the voltages» Vo and Vg were held at ground | 


"potential. However, the voltage Vp was a functionof time and forthe particular. 


j “ 


_ problem studied varied sinusoidally. This voltage was generated by the low. 
es frequency function generator and, as before, the output aes of all ampli-— 


om Case III.—For this problem the boundary at is moving with some 
scribed velocity Up = Uo folt) and the pressure —=0 -There-— 


pony the voltage, Vo, was” made to vary in the same manner as fo (t) with Vp 
0. For the particular problem considered the voltage Vo varied sinusoidally. 
The difficulty encountered inthis ‘case wasthat y y infinity andthere- 
_forethe nodes could not be equally spaced. ee, an inspection Eq. 17 in- 
- dicated that if 7 # 1 then the Up + 1, Um - 1, Um terms must be modified by 
factors that are functions of n. In the computer circuit this ‘corresponds to 


_ modifying the input vc voltages by these factors. This can be readily accomp- | 


— 
x 
, — 
m ¢ ¢ ¢ | 
| 
— 
4 
FIG, 4.-BASIC COMPUTER CIRCUIT 
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lished / through the us use of the coefficient potentiometers. As previously men-_ 
tioned, it was desirable to make 7 very large. For the particular problem | 7 


studied, the actual value of n did not appear to be too critical for this last in- 
terval 


THEORETICAL SOLUTION 
- gain some insight into the accuracy of the computer § solutions, ‘it was 
thought advisable to compare > these solutions with the theoretical solutions. 
_ The basic problem » to obtain a solution to Eq. 1 subject to certain boundary 
and initial conditions. The general solutions are presented here, and the speci- 
fic problems of interest are treated as special cases. For Cases Iand II, in 


which the boundaries are fixed, the solution can be formulated in the following 


was assumed that. & pressure gradient ‘was a per iodic function of time Boe 


oly) cos nwt + ) 


It then on 


+ + Bn Caly)] cos not 


i 
[By Snly) - An sinnw 


in n which is the absolute viscosit and 
= q h + sin y) sinh q(bty) (18) 


cos — y) cosh nat? y) + cos at» y) cosh eins y) cos 2qb - cosh qb 


2qb + cosh 


6 “Critical Numbers for Steady and Pulsating Flow,” L. J. Kastner 
= Ss. H. Shih, Engineering, Vol. 172, 1951, p. 389. 
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3 For the special case ee in the present investigation 


+ + Ay cos wt ‘Bh = eee 


w 


2qb 


The transient solution is obtained from the equation 


the initial condition uly,o) = 0 
‘The sc solution g given _— the sum of Eqs. 17 and 26 thus represents the general 
"solution tothe problem for which the pressure gradient is an arbitrary periodic 
_ functionof time; and the sum of Eqs, 22 and 26 gives the solution for the special 
— in which the pressure gradient consists of a constant component that fluc- 
JA computer solution was first obtained fort the transient flow developed when 
_ This type of 
flow is commonly called “starting flow” ‘and is s referred to herein as Case L 


~The theoretical : solution for this case is obtained from lead . 26 in which the coef- 


| 
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Sty) = + cosh 2qb © «°° 
os 2qb- cosh 2qb 
| 
by separation of variables, This solution is 
| 
— 
4 
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6 ANALOG LOG COMPUTERS 


ficients are obtained from Eq. 27 
al 


Ao 


A comparison of the resuks. obtained by the solution with those 
- obtained from the computer is shown in the dimensionless plot of Fig. 5 for > 
three representative times. agreement between the two solutions appears 
‘to be satisfactory, mits 
_ As noted previously for the problem in which the pressure gradient is sa 
function of time the general solution is given by the sum of Eqs, 17 and 26, 
Eq. 17 indicates that the “steady” part of the solution consists of an n oscillat- : 
ing component of the velocity plus the characteristic 
solution olution was obt for ‘Case Iw with 
Acomparisonof the computer solution with the theoretical solution is shown 
inthe dimensionless plot of Fig. 6. Only the results for the steady oscillating 7 
part of the solutions are shown, However, it should be emphasized that — 

_ transient part of the solution can be readily obtained from the computer, The — 
parameter qb was arbitrarily selected, and the velocity, profile at the 
when the velocity at the midpoint nenmoen two plates isa maximum is shown 

in Fig. 6. As in the previous case, the two solutions agree closely. 

a The theoretical solution for Case III can be formulated in a manner similar 

- to that of the first two examples, For this case the fluid extends to infinity in 
the positive y-direction and is bounded by an infinite plate having a prescribed a 
_ periodic motion in the x-direction, In general this motion can be specified as 


008 nw 


= on (y) ¢ cos wt (y) s: 


(32) 
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| - 
= 
_———_ . The other boundary condition is u(#, t) = 0, and the steady state solution can &g 
| 
= 
COs 5 y + Cn sin 5 y | Sin 
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© Computer solution 


FIG. _ 7.—COMPARISON OF THEORETICAL AND 


FOR CASE II. wii 


The e transient part of the can obtained ina a manner similar 
that described for Case . However, after an interval of time this part of the 


For the special case , with Dn = 0, Cn = Uo, and. a® "4 the steady state so- 


73 
(cos 


= 


A dimensionless plot of the velocity profiles for the times when the velocity 


at y =o is a maximum andwhen equal to zero is shown in Fig. 7. The parame- 


ter r is the wave length and is equal to 2 ™/ 5 The computer s solu ution is also 


‘shown, The agreement between the two on... is ; reasonably gx good, but not: as 


_ Although the specific examples considered inthis investigation are relatively , 
- ‘simple, the same technique « can be used for more involved problems, that is, 


iia 
— 
— O02 O04 4206 10 
— 2 
4 
(34) 
7 
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fo or more general types of pressure gradient fluctuations or boundary motions. 


As far as the computer solution is concerned, the only requirement is that a 
voltage to simulate these variations can be generated, Of course, as noted 


from the general solutions presented, theoretical solutions can also be obtained 
for more involved problems. However, the major advantages of the computer 
solutions lie inthe ease with which the various problems canbe solved and the : 
_ parameters changed, ar and in| the complete picture of the solution which is ob- 
tained, that is,a ,a continuous recordof the e velocities, including both the transient 
and steady state components, The results of this study (as illustrated by Figs. : 
_ 5,6,and 7) indicate that a reasonable degree of accuracy can be obtained with _ 
a relatively simple computer circuit and by using a small number of intervals - 
in setting upthe finite difference approximation, As the problem was set up on a 
the computer, one amplifier per node was required, However, for problems 
& which symmetry exists, there is really only a need to solve the problem for — 


a portion of the field.3 For example in Case I it is recognized that, due to the 
of the the velocity distribution is symmetrical with respect 


‘i , technique discussed herein can be readily extended to uniform two- 
=_—" laminar flow problems. _ The basic differential equation for this 


_ difference expression. For this type of problem it is important to take ad- 
vantage of symmetry so that the number of amplifiers required ¢ can be held at 
=. a reasonable figure. This procedure has been successfully applied to two- 
dimensional heat transfer problems.? Although theoretical solutions to Eq. 35 
are available for various boundary geometries, they are rather involved, and | 


mentioned previously, the s simplicity and flexibility of the computer solu-— 


tions make this approach particularly attractive, 


before the left. can approximated by ar an n appropriate 
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APPENDIX.— _NOTATION. 


The following symbols are adopted for use in this paper: 
= Coefficients of ‘Fourier series used to describe pressure gradient; 


4 
— 
4 — 
a 
tm iim 
pore 
ed by the Iowa State 
a 


- 


= onaiteanadtne of Fourier series used to describe boundary motion for : 


= CO ynstant of proportionality between velocity and voltage; 
constant t of hongenmenine between actual time and computer time; = 


= - frequency parameter q= 


= resistance; 


velocity; 
= of ‘boundary; 


= steady state am 
transient component of 


voltage; 
= variable input 


= ltiplier for t voltage; 
= coefficients of Fourier ‘seri es us ed. to de describe st eady state velocity ; 


for CaseslandI; 


= fre uenc arameter 


= multiplier for increments | Ay; - 


= wave length; _ 
= kinematic viscosity; 


= ‘computer | time; and ‘a 
= angular fr frequency. “J — 


— 
Pn * 
n@ 
© 
q 
— 
Py coefficients of Fourier series used to describe steady state velocity 


= 


4 


ENGINEERING MECHANICS DIVISION 


a of the American Society of Civil a 


BENDING ANDS STRETCHING | OF LAMINATED AEOLOTROPIC PLATES@ 


Yehuda ‘Stravsky! pad 


‘subject to plane forces and loading. It i is shown that a 

phenomenon occurs, ‘between in- -plane- forces and transverse bending, in the 
a 7 differential equations as wellas inthe boundary conditions. Some special types 

of heterogeneity are indicated for which the coupling of stretching and bending = 

- disappears. The general ‘results of the paper are illustrated by means of ex-- 


aq In the following, a linear theory is established for ‘the bending a and stretch- 


ing of a thin, omen, aeolotropic, multi- layer plate. ae = 


plates a coupling phenomenon. occurs between in-plane stretching and tra trans- 


~ Note. —Discussion open until J May 1, 1962. To contend the closing date one month, a © 

written request must be filed with the Executive Secretary, ASCE. This paper is ae 

_ of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the 
_American Society of Civil Engineers, Vol. 87, No. EM 6, December, 1961. oe 
_ @ Revised and extended version of paper presented at the October, 1960, ASCE An- 

_ nual Meeting at Boston, Mass., based on a part of athesis presented to the iad: 
Inst. of Tech., Cambridge, Mass., in September, 1959, in partial fulfilment of the ne re- 

1 Senior Lecturer, Dept. of Mechanics, Israel Inst. of Tech., Haifa, Israel; formerly, — 
Prof. of Structural Engrg., M. I. T., Cambridge, Mass. 


2 “Bending and Stretching of Certain Types of Heterogeneous Aeolotropic Elastic 
‘Plates, ‘Relssner and Journal of Applied Mechanics, No. 61, APM, 
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similar coupling was previously observed? for an asymmetrically lami- 
usted orthotropic shell, At the time the present paper was written a relevant — 
paper appeared,4 in which the coupled differential equations were derived for 
a multi- layer isotropic plate } ING non-linear strain-displacement rela- 
tions. It was concluded that “in the small deflection theory the in-plane forces | 
_ and transverse deflection are incoupled. . . solutions may be obtained directly 
"from the single-layer plate by transposition” = | 
_ The linear theories presented elsewhere2, 3,4 are obtained as special cases 
a = the general theory developed herein, and can be applied to the analysis ofa _ 
_ plate laminated of aeolotropic, orthotropic, and isotropic layers. = 
ae It becomes clear that even for the small deflection theory of K. S. — 
and Ss. B. Dong, 4r no simple “transposition of single layer plate theory” is pos-— 
_ sible, because stretching and bending are still coupled in the boundary condi- 


tions, although thi the « differential equations are are discoupled. 


= OF THE PROBLEM 
Consider a a composed of homogeneous aeolotropic s sheets, 
conveniently, t the bottom sheet is coe ngige as “layer 1,” the top sheet as “layer 
and the origin of the coordinates is 3 chosen to coincide with an arbitrary 


Taking the positive direction of the thickness coordinate, Z, upward t the ge- 
_ ometry of layer m, for example, is defined i by its bottom plane, z = hm-1, and 
z= hm, with thickness (hm - hm- m-1). The total thickness of the 


a Defining stress resultants, stress couples, reference surface strains (at 
= 0), and bending curvatures, as is usual in plate theory, the following equi- 7 
‘Hbelum equations and strain- displacement relations ‘of linear plate theory, 
based on the Euler-Bernoulli h othesis are ossible: 


= 


@ 


“The Calculation of Laminated Shells,” by S. A. Ambartsumtan, wv. 
Akad. Nauk ArmSSR, Ser. Fiz. - Matem., Estestv Tekhn. Nauk 6, 3, 15- 35, 1953. oe: 
*Blastic Bending of Layered Plates,” by K. S. winter and S. B. nd 8. B. Dong, Proceedings, 


ASCE, Vol. 85, No. EM 4, 1959, 1. 
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| ae 
Gy 
| | 
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i 
= 0 , M +M - Q = 0. 

e9+zK =€9 +zK .... (4) 
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AEOLOTROPIC PLATES» 


7 _¢ The components of strain, according to S, G. Lekhnitsky,5 define compo- 
nents of stress through stress strain relations that can be written, for example, 


for layer m, using a convenient notation for the elastic moduli in the form: © 


ee 
the eunmination E. Reissner and the writer,2 it is as 


that E™ varies arbitrarily from layer to layer but follows the symmetry rela- 
tions EM = EM (i,j = x,y,s). a 
In order to obtain a complete system of plate Eq. 7 is introduced, | 
in view of and 6, the N, 


The elastic | areas, elastic statical and the moments of 


az, 


35 “Bending ofa a a Nonhomogeneous Anisotropic Thin. Plate of Symmetrical C Construc- 


tion,” by S. G. Lekhnitsky, Journal of Mathematics and Mochantes, P} 


lc 

mm] = |em pm pm | | 
computed with regard to the reference plane,z = 0. 
hee of plate stress-strain relations is obtained as follows: a 4 
x) | “xy “xs xx xy x 
" 

: 
gg 
= 


‘It is clear that whenever E.. are functions Z, no summation will 


be necessary in Eq. 9, and the integral limits willbe z = Oandz=h 


8 also_ be written w with | of suitable matric ces 
‘the form 


Eqs. 1, 2 and 3 together with Eq. 8 constit rute mot eleven equations 


for eleven unknowns N_, Qs Myr 


"function F, defined 


is introduced ‘next and assumes that the body forces p Pxs Py in Eqs. 1 are re- 

stricted so as to bederivable from a potential, 2, as a = m 

ed ae the same method that appears: in section 3 of the e work of Reissner 
and the -writer,2 3 and are reduced to two ‘simultaneous fourth order 


oo 


2A 


2 


4, 
L= pe —— +2 (D +2D ) 


By y 


It is) noted that although A* and D* are symmetric matrices, the iin need 
true for B”*, Eqs. 12 17 are remarkable for two specific prop- 


r ing that is governed by the L3-operator, the structure of which is con- _ 
nected with the B-matrix inthe stress strain relations in Eq. 8. Ingen- _ 
-eralization « of the type of heterogeneous plates analyzed by Reissner and 
writer,l all the elements inthe B-matrix are regularly obtainedand, 
consequently, all terms appear in the L3-operator. This means that 
tha case, evenfor one dimensional problems, w and F are still coupled, — 
operator, L1, involves me elements of the matrix, D*, rather than 


— 3 
| 
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Four fp(x x,y) = 0, are associated with 


the system in Eqs. 12, and take the following form when ‘displacements or 


~<a In view of Eqs. 5 5, 6, and 8 the coupling of wand F enters, 1 not ‘only into the 
differential Eqs. 12, but into the boundary conditions of Eq.19 aswell, | 
Once the system of Eqs. 12 is solved and subjected to suitable boundary 7 
of Eq. 19, allother unknowns canbe determined directly. The stress 
resultants are given by Eqs. 11, the couples by Eqs. 4 through 6 and 8, and 
the transverse shear resultants by Eqs. « 2. The strain components are one 
by Eqs. 4 through 6, , and the stress components, are re obtained from Eqs. te a 
GENERALIZED STRESS FUNCTION 


When body forces are not considered it is ‘convenient, in some problems, 
to transform the system of two § simultaneous differential equations for | wand 
“Fa as in Eqs, 12 intoa single — on of the eighth order, in terms of a gen- 


(21) 


w, = or M =M_ 
il 
é 
introduce into Eqs. 12, noting that the sec 
while the first one reads _ 
©. 
L= L, La + Ly*.... 
AS De. + ( 24) 
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-8A* p* +2 


yx \ ss) 
At + + at ) p* = 8 A* D 
+2 Be +B* -2 2B 
=-2a* p* ( 


XS XX 


+4A* De +2 Box 

“ys at 


+ 2 Ds, 
xs \xy 


B* 
Bex) 
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sy Boy 


‘similarly, the the boundary conditions in ‘Eq. 19 can be expressed. in terms of 
so the problem is completely defined. 


wens 


CYLINDRICAL BENDING OF LONG RECTANGULAR PLATES © 
I 


nsight into some of the | consequences of the general elastic meee 
_ considered here may be obtained by determining the state of stress and defor- 7 


4 A* i (D +2 D ) 
Pe and | 
9 “XxX = xS “xS \ xy 
= UTC U y — 


Consider plates of layers that follow the stress strain relations 

7 given in Eqs. 7. Let the y-axis coincide with one of the longitudinal edges of 
the plate, and let 1 be the free span of the heterogeneous plate. . The edges x 

= 0,1 are assumed to be simply supported, and the edge x ° = lis free to move ) 
in the x-direction while the other one is fixed in place. wm rails _ | 


3 Supposing that the plate is also supported at y = +00, it is clear that the be- _ 


‘havior of the laminated structure under the load p= p(x) will be independent | 
y; the problem is, therefore, one-dimensional, 

This simplification pereets on an im immediate integration of Eq. 21, which takes 

now now the following form: 


solution of Eq. 33 is is ‘subject the following conditions: 


Integrating E 


view of of Eqs. 20, 4, +5, 6, 8, ‘and 34t the following is obtained: 


[ot -}- 


— 
ay 
Ky 56x 720 ~ Tx 720 720° 120° 
4 
| 
mM. = Bt | 
My = Byx Byy x(1 - x) (41) a 
= Bfs + Dis x).......) 
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| Making use of ‘Eqs. 4, 5, 6, 8, 36, and 37 the strain components 2 are: 


7 


ess. components, in view of - , 43, 4 


* 


‘The are obtained by using Eqs. 2, 40, Ai, and 


AY 


= Bie Ai Ba] 
macroscopic equilibrium equations, which for the problem takes the f oe 


— — 
Nem (at, Byx - Ajy 
ty = (Aly Ajy Bis) 
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in terms of the transverse | shear it is found 


i=l 


Fin 
228 
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xy ‘yy xx 
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Note that in order to maintain the deformed shape in Eq. 36, a cross re- 
sultant force, Ny, will be. developed in the ae in addition to the couples Mx, 
My, Mxy and their shear resultants, Qx, Qy, that occur in the homogeneous 
-aeolotropic plate ‘subject to cylindrical bending. 

_ Eqs. 36 through 56 are also remarkable for the state of stress and strain. , 
that characterizes the behavior of the heterogeneous pute. . Transverse bend-— 

- ing is associated with in-plate deformations €, and | €xy and in plane stress a 
Ty, xy of quite a complicated dependence on the elastic ‘moduli, elastic 
elastic moments, and moments of inertia of one 

a 
UNIFORM DISTRIBUTION OF STRESS RESULTANTS AND COUPLES 
Assumed that the plate is subjected to resultants, N, and couples, M, that 
are independent of x,y, an attempt is made to solve the system of Eqs. 12 | 
_ the rectangular plate bounded by x = 0, = and wae 0, sb — to the following — 


boundary conditions: 


-= an 


m,R =0 


Because all prescribed fo forces a and are constant, itis established, 


view of Ea. 10, when Nx, Ney have constant values and 
there are no body forces ‘Px, Py 


In order that Mx, My, Mxy at the same time have constant values mx, = 

Mxy it is required that kx, Ky, Kxy have constant values kx, ky, kxy. -Accord- 
ingly, the deflection function, w, takes the form Ope 


> » ot 
— 
&g 
: 
— 
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EM 


in a, a 

following relations: 
a 


in in which 


Eqs. 61 57 it: is that 


“through 7 in conjunction with Eqs, 59 and 60, 
_ The expressions for the stress: components in terms of stress ne 


and couples are mentioned here: = 


Com] = (em) ([a] + + [em] + 


_ which expresses a remarkable characteristic of heterogeneous plates: gener- 7 
ally, each stress component is a function of all and congien, 


_ SOME SPECIAL CASES DERIVED F ) FROM. GENERAL THEOR' THEORY 
‘The specialization « of the general system of Eqs. 12 to the —_ all 


isindicatedheres 
I, For which bending and stretching are uncoupled 


IL For which the cross- elasticity effect, Eq. 70, ae — 


ud To determine the constants in 
(62) 
= 
nmetric matrices but b, c are, generally, not symmet- 

and, thus, th re uniquely determined. 
ctrain and ctrescc contribution can he comouleted hy using Kos 61 and 4. 
ff]... (70 

— 


An example of case I is the laminatedacolotropic plate, Let 

the plate described in section 2 be the upper half of a composite plate of total 
ma, h, made up of 2n layers, For such a plate it is convenient to take 

the: reference plane z =0 to coincide with the symmetry plane of ‘the plate, | 
“In view of Eq. 9 it is established th that ‘13 = 0. The system Eqs, 12 is uncoupled 


el 
and takes theform 
4 Ld A71) 


eevee eee 


in which 


ee w , and F goes. This was first shown by Lekhnitsky5 for the bending ie 
 shownin Eq. 
‘However, it is noted that the effect in 70 is still q 
but is simplified to | 


= [am or [rm] =) 


76), 


and D is given by Eq. 74. 
As an example of case II consider the multi- layer aeolotropic plate, for -_ 
which all the elastic moduli E in the manner with z . Thus 


ail 


4 2 é 
: | : _ The difference between Eqs. 71 and 72 and the classical equations of homo- a 
; | 7 _ geneous aeolotropic plates lies only in the value of the constants in the opera- ig 
a | tors Lj and L2. Known results of homogeneous plate theory can be fully uti- i i 
— 


that 


in which I is the unit matrix. 

_ Because B* is a diagonal matrix with equal elements, the coupling operator 7 
vanishes identically inviewofEq.15, | 
quantity, given by Eq. 80 is identified as the. distance of sur- 
_ The existance of a single neutral plane, at which allstrainc components ° van- ; 
ish simultaneously, causes the cross- elasticity « effect of Eq. 70 to disappear, 

= 0 at the neutral plane of the 


Eq. 61takes the form 


— 
— 
4 4 em ( > — 
— 
| Using Eqs. 82, 78, and 70 it is found that ™ oe . 
mile om [ 20] (fa — 
[rt 1] or } [a] LN] 


_AEOLOTROPIC PLATES 


are given by 


is noted ‘that, although cases I and II are mathematically identical, case 
a II is completely analogous to that of classical stretching and bending of homo- 
geneous aeolotropic plates, but case el involves cross- -stress systems, as in in Eq. 


For example, bending moments Mx causes self- stresses S Oy on 
boundaries y = 0, ly. 
_ The far reaching ¢ consequences of the presence of the uncoupled cross- 4 

- “elasticity effect (Eq. 75) in case I show up ina plate theory that considers 
_ transverse shear and normal stress deformations. However, an examination — 
of this problem is beyond the scope of this paper. a ~— 


To obtain a insight into the resus of uniform distributions. of 


=1, -1i in Eq. 


[- 2.7822" 
y2 + 6, 6.737039 x 


2, 782271 


“On the Theo Heterogeneous Anisotropic by Stavsky, thesis pre- 
‘sented to the Massachusetts Inst. of Tech. ., in Cambridge, Mass., in September, 3959, é 
in partial fulfilment of the pe for the degree « of Doctor of Philosophy. 
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Ana logous results for the case of ai uniform m bending soment, Me an are: 
- 0.176025 408201 x x 


+360902 0.927424 3. 


-1,978923 | 


22.3 344940, 


A graphical representation of these results is given in 1 Figs. 1 and 2, 


2.— R = = 1, m = 8 in Eq. (88. 
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AEOLOTROPIC PLATES 
CASE 
I 
Nx 
_FIG, 1.—STRAIN AND STRESS COMPONENTS IN A HETEROGENEOUS 
AEOLOTROPIC PLATE, WITH THE SAME Z VARIATION OF 
SUBJECT TO UNIFORM AXIAL TENSION 
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4 ‘FIG. 2. —STRAIN AND STRESS COMPONENTS IN A HETEROGENEOUS AEOLO- 
TROPIC PLATE WITH THE SAME E VARIATION OF OF Eyy, SUBJECT 
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Ag results for case 2 is given in 3 and 


cross- elasticity effect of Eq. 70 vanishes for the type of hetero- 


geneity indicated by Eqs. 88 and 89, only 7x . stresses are developed of quite 
__It seems striking how effectively one can control the elastic stress dist 


- bution ” the type of the heterogeneity of the plate and, thus, the material is” 


> 


> 


ZO) 


GENEOUS AEOLOTROPIC PLATE WITH THE SAME Z VARIATION OF E 
TO UNIFORM AXIAL TENSION Nx 


“utilized much more effectively in n bending than inthe conventional ¢ elastic homo- 


geneous system, 
So far, a special class of heterogeneity has been dealt with; namely, ; all 


- elastic moduli have the same variation throughout the plate thickness, A s some- 
what more general case is now considered as follows: : - 


ENTS INA SYMMETRICALLY HETERO 


0.2 (1+ + 0.6 cos Fz) 0,002(1 + 100 cos 
0.002(1 + 100 0s fz) 0 0.002(1 +100 cos 1 


— 
— 
= 
g || 

) |. (102) 


= 


0.688 


“FIG. 4.— _STRAIN AND STRESS COMPONENTS INA HETEROGENE- 
_ QUS AEOLOTROPIC PLATE WITH THE SAME Z VARIATION OF Ejy, SUB- 
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ogous results for the case a bending n moment, Mx, a are: 
‘err 
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"The following ‘results were obtained for a uniformly axial force applied at 
_ the middle plane of the plate (2 ~ 
103) 
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FIG, 5.—STRAIN AND STRESS COMPONENTS IN A NON-SYMMETRICAL HETERO- 
— | = GENEOUS AEOLOTROPIC PLATE SUBJECT TO UNIFORM AXIAL TENSION 
— 
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of the for Case 3 is given in Figs. 5 and 


It is noted that the strain and stress distribution deviate considerably, quan- 
titatively and qualitatively, from the corresponding results for homogeneous — 


_ The occurance of cross-stresses is clearly examplified and it seems that 


it has to be considered in the design of laminated structures, = Ae | 
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IN FLUID DYNAMICS". 
By J. W. Delleur,! M. ASCE, and A. A. 


“sented. The case in which the boundary conditions specify the velocity on | 
part of the boundary and the stresses on the remainder of the boundary is con-_ 
"sidered. The method is particularly well ‘suited for machine computations. 
The variational principle is applied to an approximate computation of the ve- = 
locity distribution for laminar flow in Square and rectangular ducts. Approxi- 


mate and exact solutions are compared. ; 


‘Warlational procedures may be particularly to obtain approximate 
solutions to a wide variety of problems. Variational techniques make use of | 
certain broad minimum principles. For example, in structural engineering 
and in elasticity the concepts of least work, minimum potential energy, mini-— 
_ mum s strain, and complementary energies are commonly us used. es Similar broad 
minimum principles applicable to fluid-mechanics problems a are not presented 
in standard textbooks. The purpose herein is to state a variational principle 
that is applicable | toa class of problems in fluid mechanics, and to Tofustete 


its usefulness in solving problems by use of digital computers. _ = 


= Note.—Discussion open until May 1, 1962. To extend the closing date one month, a 
- “written request must be filed with the’ Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings | of the — 
American Society of Civil Engineers, Vol. 87, No. EM6, December, 1961. ome a 

* Presented at the October, 1960 ASCE Convention in Boston, “Mass. 
Assoc. Prof. of Hydr. Engrg., Purdue Univ., , Lafayette, Ind. 
Research Asst., School of Civ Civ. ‘Engrg., Purdue Univ., , Lafayette, will 


4 
| 
— 
q 
Ly 
4 | 
— 


58 

ey ‘The basic variational principle of particle dynamics was known by L Lagrange. 
; - The Lagrangean equations of motion, Bernoulli theorem, and other equations — 


of hydrodynamics in the Lagrangean system can be derived from variational 


principles. 3 However, engineering problems are usually formulated in the 
_ Eulerian system. The variational principle stated herein is the Eulerian sys- 
in which inertia forces are negligible. 
- Important features of this pr principle: are ‘the two boundary conditions considered: © 
The case where velocity is specified along the whole boundary or the 
stresses are zero along the entire boundary; and 
The case where the boundary conditions specify the velocity on one part 
of the boundary and the stresses on the remainder of the boundary. — ites Ae 
The principle was first stated by Helmholtz4, and the discussion of the gen-— 


— point of view was first given by C. &B. Millikan, a 


eles equations \ will be reviewed before the principle is stated, proven, and 


_ Notation.—The letter ‘symbols adopted for use in this paper are defined 
ohare they first appear, in the illustrations or in the text, and are arranged . 

alphabetically, for convenience of in the 


VISCOUS- -FLOW 


problem in flow may ay be formulat formulated in terms | of the differential 
equations describing the phenomena and the boundary conditions of the prob- 


lem. The governing differential equations for fluids? are 


3 “Variation Principles of Hydrodynamics,” by c. _ Eckart, The Physics | of Fluids, 
“Zur Theorie der Stationéren Stréme in Reibenden Flissigkeiten,” by 
Helmholtz Verhandlungen des naturhistor isch—medizinischen ‘Vereins zu 


_ “Hydrodynamics, ” by Horace Lamb, Dover Publications, New York, N. Y., 1945, 


“On the Steady ) Motion of Viscous, Incompressible Fluids; with Particular Reference 


- to a Variational Principle,” by C. B. Millikan, EnSoseghion! Magazine, Vol. 7, No. 44, | 


“Modern Development ‘in Fluid Dynamics,” by s. Goldstein, oxford Univ. 
Oxford, England, 1938, Chap. 3. 


: = 
— 
tthe Naviier-Stokes equations of motion 
which ‘iss velocity, t denotes time, p is density, p represents pressure, the 
indices i, j, and k denote the x, y, and z components, and summation should 


VARIATIONAL METHODS 
made over the three values of a repeated — occurring in the expression for 


uj 


in _— bij is the Kronecker delta, equal ‘to unity when i=j and equal to zero 
= other conditions, and is absolute The rate Of of 


oF VARIATIONS . AND EULER- LAGRANGE EQUATIONS 
‘The calculus of variation makes use of broad minimum principles that char- 
acterize a class of phenomena. It is then concerned with the extremization 


_ tities on which J depends. For example, it is desired to find what plane’ curve 
-connecting two points (x1, y1)a and nd (x2, y2)t has the shortest length. The length 


the curve is given by — 


jp. 


which pe pass through th the given end potas, | the shortest one 


must be prot thus reduces to finding the function that 


makes the ‘integral a minimum. Therefore, there are two features to be 
preven 


2. The condition in addition to the boundary conditions bat 


such as y(x) must satisfy in order to make Ja minimum. 
mathematical condition is first the physical principle will 


6 


FG, dx . 
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FIG. 


~ 


the function y(x) boundary conditions and the Lagrange 


7 In fluid ‘mechanics the velocity may be a function oft two or three eindependent — 
| a variables. The previous theorem can be extended to several dependent and 1 in-— 
u, Ux, uy 


satisfy the Euler- differential equation 
= a oF oF \ 
in addition to the boundary conditions. 
“A The o object of the variational principle is to state a minimum principle | that i 
characterizes laminar flows. _ This makes possible the construction of of inte- 
-grals that attain their minimum when the flow is steadily uniform. The ae 
: lem of determination of velocities reduces to standard methods of the calculus > 
of variations. Once velocities are known, pressures and other quantities may > 
be calculated by the usual methods of fluid 
The variational principle is as follows: The motion mn of an incompressible 
4 fluid that satisfies the equation of continuity, , the equations of motion, and andthe 


attains a minimum for uniform flow. in instead, the ‘boundary condi- 
- tions specify the velocity on one part Sj of the boundary, and the stresses nd 
the remaining part So of the then the expression isa min- 


in which I is the excess of the dissipation integral over twice the rate at which ; 
work is being done by the specified ‘surface tractions oi (Fig. 2). of Hence, for 


8 “Calculus of Variations with | Applications to Physics and Engineer ing,” by R. Wein- 
stock, McGraw-Hill Book Co., New York, N. Y., 1952. 
“Methods of by R. Courant and D. Hilbert, Interscience 


Publishers, Inc., New v York, +. , 1943. 


oe 
— 
4 
ig 
4 
ary is such that the dissipation integral _ 
— 
| 


work is being done on the boundaries over the total rate of change in kinetic — 
plus potential plus dissipation of of energy ' within the e boundaries 

7a variational principle may ‘be: justified by replacing ‘the integrand Fo 
the variational integral in the Euler-Lagrange equations corresponding to _ 
_ variational principle stated previously, and by showing that the result obtained | 
is equivalent to the equations of continuity and of motion describing the prob- 
lem. In this fashion it can be formally proved that minimizing the integrals of 
“ertia terms neglected: ‘this fact is used in 5 the solution | of problems. If the 
closed-form solution of the Gorey “ motion and continuity fora given prob- 


tial equation does not need to be solved rigorovaly. Approximate methods 


available, among which are the Rayleigh-Ritz and the Galerkin methods, 8,9 the | 
former being illustrated in the examples. + 


Mia: The variational principle stated previously may also be verified by taking 7 
the variation of J orI and investigating the conditions under which they vanish. 
Taking the the variation of the dissipation integral (Eq. 1 13), one obtains 


ey +8 eij ‘ey oy) dv 6 ejj dv 


because of the oo - the stress and strain tensors. . Making us use of E 


6 uj 


steady, incompressible laminar flow, the velocity distribution is to be chosen 
hat the boundary conditions are satisfied and excess of therate at which fix 
— 
| ; 
5) 
+ 
Assuming steady uniform flow, there 1s 
o and the equations of equilibrium10 app — 
10 «Theory of Elasticity,” by S. Timoshenko and J. N. Goodier, McGraw-Hill Book 
mm IIc., New Yk, N. Y., 1951, Article 


ae TIONAL METHODS 


theorem 
J=2., @ al 6 uj dS =2 


in which ai is the direction cosine of the outward normal to aol surface ele- 


ment. The surface integral in Eq. 20 will vanish it 


or if 
= 


7 on Ss. If either condition in Eqs. 21 “e satisfied 


and the integral of Eq. 13 is an extremum. 
7. If the boundary conditions specify the velocity on one part a ‘the | boundary 
7) Sj and the traction on the remainder of the boundary = the conditions a Po 


But ut because th the e boundary | conditions specify t the daeede 2 


we S15 uy = 


ever, Eq. 20 still holds and may be 


_which shows that the integral of Eq. 14 is an extremum for mixed boundarycon- _ 


ditions. 


— 
| 
is equivalent to st ms, the term— 
— 
— 
4 
= 
J % ey dv-2 J olujdS} =O ........ (26) 
: 


he foregoing pro 
I and J aes to a 


= 


distribution ‘for laminar flow in square and rectangular ob- 
tained by the variational technique is compared tothe exact solution. The same | 
technique has been applied by D. F. Hays}! to the problem of hydrodynamic — 
i Circular Pipe.—Consider a horizontal circular | pipe of ‘radius To; and let 
u be the velocity in the z-direction parallel t to the axis of the pipe. The Navier- © 
_ Stokes equation in cylindrical coordinates, 
inertia terms is (Fig. 3): 


in the z- direction, neglecting t al 


: Because of axial nines, the third term | in in Eq. 28 vanishes, and 


boundary condition is is 


‘= r,, and the pressure gradient as on walt a as s the viscosity are assumed to be 
given; thatis,C*is known, 
_ An approximate solution of the differential equation (Eq. 29) along with the 
; 7 condition of Eq. 30 can be obtained by use of the entities steps: _ 


La). 
4 Zz bs Because the boundary conditions specify the v velocity u=0 at the wall of 


rae the head loss per unit a along the conduit is given — ie 


the problem falls into the category of of mixed boundary conditions, and 


of Eq. 14 must be utilized. xy the problem at hand the energy dissipation per 

unit time per unit volume is12 
(32) 


i “A Variational Approach to Lubrication Problems and the Solution of the Finite 


Journal Bearing,” by D. F. Hays, Journal of Basic Engineering, ASME, Vol. 81, March, . 

_ 12 “Fluid Mechanics,” by | v. L. Streeter, McGraw-Hill Book Co., Inc., New York, N. 


of shows that I and J are extrema. Helmholtz4 has shown 
that unique and stable minimum. 
the: Subsequent material the velocity distribution for laminar flow incir- 
= 
31) 
— 
> 
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| 


\ 
\ 
— 


in 1 which | So is the cross-sectional area at the en 
pressure forces ee an element of unit length in the z-d 


ds of the element on which the — 7 
I 
> z-direction, I may 


‘in n which the subscript r denotes derivative re with respect to r. — —_ 
- _ The variational principle states that the integral I of Eq. 35 | isa mini- 
mum. The velocity u(r) that minimizes Eq. 35 must also satisfy the Euler- 
Lagrange differential equation 5 (Eq. 10) which in circular coordinates become 
Performing | the of operations indicated in Eq. with the value of F given in Eq. 
one obtains 
which after dividing by 2 u reduces to This shows that the 
‘differential equation (Eq. 29) with its boundary condition (Eq. 30) gon nell 
to finding the function u(r) that minimizes Eq. 35. | This ‘is s the fundamental 
principle on which the variation technique is based. 
_ 3. A velocity distribution that satisfies the boundary condition (Eq. 30) is 
assumed. _ This is usually a polynomial or a truncated series. _The variational 
_ technique yields: the best ‘set of coefficients for the polynomial or truncated © 
series assumed. Further discussion on the selection of the assumed form 
the distribution will be given with the next example. 


og The rate at which work is being done by the specified pressure forces per unit | ae 
: The integrand of the foregoing variational i integral is an ie 
— 
fi» 
“4 
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‘This function ‘obviously satisfies the ‘condition of Eq. 30. 
Ai - ALG 
= 


* r (Ag + 


which may be integrated to. the variational integral 


ry, 


olving 


this is is ‘the well- known parabolic velocity ‘distribution for laminar flow in pipes. 


g 
A, R®) . . (39). 
— 
5. Differ (Eq. 38a) and equating to zero lution will yield the 
sumed function q the unknowns. The so 
— ions with the coefficients as the un _— 

— 


Ps the eeu example is trivial, and the variational technique is 

- not the simplest method for obtaining the desired result, it illustrates the se- 
quence of steps to follow. The method is actually a more powerful t tool when 
the exact solution of the differential equation cannot be obtained or is extreme- 
ly laborious to obtain. The variational techniques will yield the polynomial of | 
best approximation of the formassumed,. 
_ Rectangular Conduit. —Consider a rectangular duct of sides 2 a and 2 b(Fig. | 

4), ). It is desired to obtain an approximate calculation of the velocity distribu- —— 
for laminar flow of a fluid of absolute viscosity u. 
The differential | equation governing the viscous flow in a rectangular duct, a 


- obtained by simplifying aa Navier- Stokes equations (Eq. 10) iS 


gral of Eq. 14 dissipation function (Eq. §) to 


Thee dissipation integral (Eq. 14) becomes 


z-direction may | be written as 


ft which the the ‘subscripts indicate partial derivatives. 


= 
> SOlution Of Ine Drovlem Of step: 
(45) 
“ 
— 
J +ul+ 2 Ce ul) dx dy......... (49) 
i 
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~ 


J (60) 
F in the Euler- Lagrange equation (Eq. 12), , Eq. 42d r may be obtained. 

- Once more this shows that solving the differential equation (Eq. 43) and the 
—, conditions (Eq. 44) for u(x, vy) is equivalent to finding the function a 
u(x, y) that minimizesEq. 49.0 

Velocity distribution that satisfies the boundary conditions is assumed. 
_ Depending on the problem at hand the function u(x, y) may be assumed to be- 


‘7 algebraic series, or a trigonometric series. The Ritz method chooses u as 
a series of linear terms each of which satisfies the boundary conditions. That 
© 


ug + ay Uy + ‘ 
_ For the the rectangular section, possible trial functions would be 


eaeral, if the | conduit is bound by ¢ a convex polygon, of which the sides have 


+ + By y (C1) (Ag x+ 


oe For: the easinies to be — herein a velocity distribution of the t type > of 


Eq. 52 has been 1 selected: 


(x2 a2) (2 p2) 
— 4 The: function F of Eq. 50 and the integral I of Eq. 49 are e then determined. 
Differentiating I with respect to each of the coefficients contained i in Eq. 57 
and equating to zero one obtains aad ‘simultaneous equations 


te 


ree 


| 
| 
the 
Anx + Bny + Cn =O = 1 
2 


Mg A. * 
(59d) 


B+C cr2) (61c) 


ce) 
‘qs. 58 written in detail are 
ory 
— 
oo 


Mg => ray 


2) 


which r= =b/a and and in which the ina of integration wie in — - 


g the integration of Eq. 49 are: 
De= 0. 0.01015873 = 0. 


= 0.05587300 0. 00338625 
000846561 = 0.00266062 


- tron Digital Computer ben several values of re The ‘obtained values of the « co- 
The velocity distribution 


Ay ,A are presented in Table 1. 


-0.29190673 =) 05923733  -0.04665263 -0.08542354 
11216367 | —-0,00688249 | -0.01755341 | -0.00617324 


-0.02974420 ~0,00058096 | 00290610 -0,00023508 
| 0.00710333 -0.00004145 -0,00027111 -0 00000536 
-0,00309772 | -0.00000845 | -0.00005860 | -0.00000051 


and AS for the de- 


be obtained by substituting the values A 
sired depth- to- width ratio 


which was from Eq. 57 by introducing the coefficients of Eas, 
The velocity may be computed 
= 


| 

— 
4 
| 


exact 


exact 
method 


variational 
ike method | 


1.00000000 
0.997 70557 
0. 99100148 


2 07515550 


}2,09719040 1.00000000 
2.07101770 |2.09237840 | -0.99800605_ 
205857290 | 2.07831860 .99200899 


—-2.03772630 2,05570450 
00832010 | 2.02470940 
11, 98464360 


0. 94939055 
0. 92662016 


0. 89931395 
-0.83027230 0.82129386 

0.78804982 | 0.77809407 
0.74031746 | 0.72945862 

0,68675561 | 0.67469802 

—0.62701422 | 0. 61361677 4 

-0,56071281 0.54655306 

0.48744039° 0.47376169_ 

0. 82773494 0.40675562 | 0.39468755 

0.64584733 31818667 0.30795837 

0.22123129 | 0.21243844 

0.22803567 0.10873389 

0.85574998x10-7; 0,00000000 0. 40804592x10- 


.97013300 
.92288090 93449860 
.87381 


.866 21620 7 0. "39948998 
.79972850 | 1.80296170 


72294410 
q 


= 


72240950 
.52981360 
41497020 
-28687110 
14622580 
99356904 


2. 
1 
1 
1 
1 
1 

1,63532580 
1 
0. 


42512460 
30115200 3 


16356620 
01151460 
84408114 
66028682 

0,45908934 
0,23938325 
-0.00000000 


* 4 


+. 


(66) 


— 
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| 
( | 
=b/a=1,0, Ro = y/b = 0.0, valu 
C8 A +2 A, +a A, ré +o A, re) .... (4) 
the velocity distribution m 
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Similarly, , when using u — 


Umax 


~ 


The velocity distributions along the horizontal axis = 0) were deter- 
_ mined by means of Eqs. 66 and 68 for several depth-to- width ratios r. _ ad- &g 
dition n for the > square conduit the velocity | distribution was aiso computed for — 
7 several values of R2. The calculated values were compared to the exact solu- | 
- tion, 3 and results are plotted in ‘Figs. 5 and 6. Some of the calculated results _ 
_ «are also summarized in Table 2. _ Figs. 5 and 6 show that for the square duct 7 
_ (r =1.0) the agreement between the approximate solution by variational meth- 
od, and the exact solution is good for both u/U and u/umax plots. With in- 


creasing values of r, increasing deviations of u/t appear, whereas S u/Umax 
remains in better. agreement. 


_ EXTENSIONS THE VARIATIONAL PRINCIPLE 
“The v: variational principle may be extended to non- -Newtonian 
case of a simple boundary condition. R. B. Bird!4 has shown that a siete 


procedure may be used for non-Newtonian fluids, with the same assumptions = 


of steady flow, and negligible inertia terms. _ In the case of a non- Newtonian 
incompressible fluid the value of F in Eq. ‘13 i is given by Be 


for for several, models are given by Bird, W. Stewart, 


a 


ash 
_ The variational principle stated herein is applicable to the solution of the 
Navier- -Stokes equations when the inertia terms may be neglected. The meth-_ 


13 “Fluid Dynamics and Heat Transfer,” by J. G. Knudsen and D. L. Katz, MeGraw- 
Book Co., Inc., New York, N. Y. , 1958, Article MeGraw= 
_ 14 “New Variational Principle for Incompressible Non-Newtonian Flow,” ” by R.B B. ge J 
Bird, The Physics of Fluids, Vol. 3, No. 4, July-August, 1960, p. 539. 
15 “Transport Phenomena,” by R. B. Bird, W. E. Stewart, and E. L. Lightfoot, . aoe 
“wiley | and Sons, Inc., New York, N. Y.,1960,p.103, 
_— 
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tions governing the problem cannot be Obtained. sample calculations how 
. that the method is particularly well suited for ‘machine calculation, and that — 
_ reliable solutions can be obtained. The method gives the polynomial of best — 


of ‘the form for the velocity distribution. 


al 
‘The writers are to Paul F. Chenea, of the School of 
ical Engineering at Purdue University, - Lafayette, Ind., whose ideas have “a 
tributed substantially to this paper. writers are also grateful to D 


valuable suggestions. calculations for the velocity distribution in square 
and rectangular ducts were made on the Datatron — at Purdue 
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= half depth of of nth 


variational integral; 
ion, 


subscripts in indicial notation; 


Myg = 
prensa intensity; 
aes: 


di l dinat 
imensio ss co 


coordinate in pipe 2 depth- -to- -width r ratio 
in rectangular duct; 


| 
— 
— 
following symbols have been adopted for use inthis paper: 
My,Mp, . 
3 


= 
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= pipe radius; 


velocity; 


average velocity; 

volume; 

= coordinates; 

cosine; 

= = variation; 


= parameter; 


= artiteary abl 


= tangential coordinate; 
= absolute viscosity; 
ty; 


= density; 
= stress tensor, at. 7 
= surface traction; and 


= rate of dissipation mn of energy per eeu sansa time per unit of 


= Kronecker delta = 1 if i=j, =0 otherwise; 4 
> 
6 
4 
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UCKLING BEHAVIOR ABOVE. THE TANGENT MODULUS LOAD 
Bruce G. Johnston, 1F, ASCE 


SYNOPSIS 


- ‘The inelastic buckling behavior o: of a concentrically loaded strut having a re- 


duced rectangular cross section as mid- ‘section is described in detail. The | 
buckling model is similar to that originally 1 used by F.R. Shanley.2 .2 Inthe Shan- 
ley model two localized points of area were assumed, whereas in the present 
modela solid rectangular cross section is introduced which permits a detailed 
exploration of stress distribution across the section. At various load levels 
‘the stress distributions that were described intuitively by Shanley, in his orig- _— 
inal paper, are determined quantitatively. The behavior of struts held so as to 

remain above the tangent modulus is also studied, as wellas other 


inthe shape of the stress strain. curve is demonstrated in in connection with ivar- 
ious aspects of inelastic | column behavior. = 


INTRODUCTION 


ii The results presented | herein pertain toa series 5 of simulated ¢ experiments 


: on ‘structural aluminum alloy struts of various lengths. Stress distribution | 
= the section, load deflection curves, and other information are deter- 


_ Note.—Discussion open until May 1, 1962. To extend the closing date one month,a __ 
written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the 
smanqenen Society of Civil Engineers, Vol. 87, No. EM 6, December, 1961. shepunte nella 


1 prof. of Structural Engr., Civ. Engrg. Dept., Univ. of Michigan, Ann Arbor, Mich. _ 
2 “Inelastic Column Theory,” by F. R. Shanley, Journal, Aero Science, Vol. 14, No. ae 


1947, p. 261, including discussion ad T. von Karman, 
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: mined by use of a computer. There are many advantages in simulated tests, 
carried out with the aid of a computer, in comparison with real tests in anac- 
tual testing machine. No machining is involved, no materials need be acquired, — = 
and there is no scatter inthe test results! Moreover, the precision of ee 
a although based on a simulated and idealized material, permits a study of de- 
tails of behavior that is not possible in ordinary laboratory tests. It would - 1 
impossible to completely duplicate the observations that be on 
basis of the simulated tests reported in this paper. aaa a | 


“JASINSKI 


VERIFICATION 


SHANLEY 
ANLEY 


1946-1947 


FIG. _—EV OLUTION O OF THE COLUMN FORMULA 


‘The maximum load (also termed “Shanley Load”) is determined quantita- 


tively for a concentrically loaded strut that starts to bend at the tangent mod-_ 
ulus load. In properly examining the work of Shanley, 2 it | is essential that the | 
development of the Euler formula and its modifications be reviewed over the 
period of 203 years between it andShanley’s work. ‘An excellent source for such - 


oe has been en provided by N. J. Hoff. 3 Fig. 1 summarizes the principal - 


3 “Buckling and Stability,” by N. J. Hoff, Journal, Royal Aeronautics Society, Vol. 


Aero No. 123, January, 1954. 
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BUCKLING BEHAVIOR 
aon in ies form. In 1744 , Euler pres presented his evaluation of the 
- average stress at which a slender axially loaded strut of constant cross sec- _ 
_ tion will develop bifurcation of equilibrium positions at a constant load. For — 
“many y: years Euler’ s formula was not generally applied toactual design because _ 
proof tests of str ctures indicated that columns frequently failed below the Eu- 
ler load. In 1889, Considére indicated why Euler’ s formula had not been more > 
_ useful to engineers. He conducteda series of thirty-two column tests and sug- 
- gested that if buckling occurred above the proportional limit the elastic mod- - 
; ulus, *E” should be replaced in the Euler formula by an “Eeff.” He correctly 7 
_ stated that this effective modulus should be somewhere between the elastic 
_ Independently of Considére, during the same year, 1889, Engesser suggested | i 
_ that column strength inthe inelastic range might be obtained by the substitution — 
- of Et in place of E in the Euler formula. This is known today as the “tangent _ 
- modulus formula” and has been adopted by Column Research Council4 as the — 
“proper basis for the establishment of working load formulas” for both ferrous 


In 1895, Jasinski suggested that there was an ‘apparent mistake in Engesser’ s 


in the inelastic range 2 should be considered, as had been done ina very general» 
“way by Considére. Engesser proceeded, within the same year, to produce a : 
7 “corrected” general formula for a “reduced modulus,” and he stated that this 


reduced modulus depended not only on “Et” and *“E*” but on the shape of the 


cross sectionas well, 
In 1910, Theodor von Karman, Hon. M. M. ASCE, derived explicit expressions 
for the “reduced modulus” for both the ae a and ‘the idealized H H-section 


VE 


. and for about 35 years subsequent to von Karman’s work, a controversy was waged 
_ over the comparative merits of the tangent modulus and double modulus column 
Sa. From the classical instability concept the double modulus tory 
was correct, because it indicated the load at t whicha perfectly straight and — 
trally loaded column could have neighboring equilibrium configurations with» 
Bo change in load. This is identical in concept to the Euler load in the elastic - 
, buckling range. However, many experimenters found that columns tested in i 
the laboratory with utmost care usually buckled at loads just slightly above the 7 
tangent modulus load. For example, very careful tests were made in the 1930’s 
the Aluminum Research Laboratories. 5 One of their conclusions was: 
test data presented hereinare inclose agreement with Engesser’s formula. . .” 
7 4 By this was meant the tangent modulus formula, even though Sageater himees himself 


; 4 “Guide to Design Criteria for Metal Compression Members,” Column 7 
Council, published by Column Research Co Council at at the Univ. Michigan, , Ann Arbor, 

“Column Strength of Various Aluminum Alloys, * by R. Templin, R. G. E. 
C. Bare 


Hartmann, and M. Holt, Aluminum Research Labs., Technical Paper re - a 
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In 1946, Shanley6 r reconciled the c controversy between the proponents of the 
_ tangent modulus and the double or reduced modulus theories. His explanation — 


"possible for acolumn to bend with increasing axial load , without 

strain reversal, it was reasonable to conclude that such bending would start at 

- the tangent modulus load. Thus, normally, for the usual stress-strain curves, 
the double modulus load never could be reached, because it is based on equi- 

_ librium configurations in the neighborhood of a . perfectly straight column. In 
a letter published jointly with the 1947 Shanley paper, von Karman2 redefined — 
the | tangent n modulus load in a way that “may be paraphrased as follows: 


“The tangent. load isthe smallest value of the axial load at which 
_ bifurcation of the equilibrium positions can occur regardless of whether 
or not the transition tothe bent position requires an increase of the axial 


Shanley showed intuitively that, after bending started, for an any ‘finite load 
- increment above the tangent modulus load, there must be, to attain an equilib- 
rium position, some ‘strain reversal. Equilibrium in a bent configuration above 
the tangent modulus load obviously is not possible if ‘everywhere across the 
— section the stress- strain relationships have been governed only by the tangent 
modulus. . This is true, even if the tangent modulus suffers no decrease any-— 
where within the cross ‘section. 


the | H- ‘section with flexibility over the full of the column. “The 
flanges were represented as in the Shanley model, , by two point concentrations — 


of area with an intervening web of zero area which could, nevertheless, trans- 
mit shear. The result is an important contribution to column theory, in that 

_ bending occurs under varying inelastic conditions all along the column, Duberg a 
and Wilder confirm the Shanley concept that bending is initiated at ‘the abeet - 


“if the behavior of a perfectly insane column is regarded as the limit- | 
-_ ing behavior of a bent column as its initial imperfection vanishes, the 
; . tangent- modulus load is the critical load of the column, -that is, is, the | load a: 
eal In 1950, Tung- Hua presented results a his inelastic 2 analysis of 
_ slightly curved column, including effects of strain reversal, ee a rec- 
-tangular section with distributed area. 
ar Notation. —Terms used in this paper are defined where they first a appear and 
are listed alphabetically, for convenience of Appendix, 


_ A sumerical procedure will be presented for the evaluation of the maximum 
(Shanley) load. Although the method is applied herein to an inelastic — 


6 “The Column Paradox,” by F. R. Shanley, Journal of Aero Science, Vol. 13, No. 5, 
a a. 7 “Column Behavior in the Plastic Stress Range,” by J. E. . Duberg, and T. W. ‘Wilder ~ 
ML, Journal of Aero ornese, ¥ Vol. 17, No. 6, June, 1950, p. 323, see also NACA T. Ne f 
8 “Inelastic Column Buckling, ” by Tung-Hua ‘Lin, » Journal of Aero Science, Vol. 17, 
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1 eomagater will be. required in such an undertaking, | and at 
ia successive e equilibrium evaluation, the column configuration by the New- 
mark9 numerical procedure will be determined. 

_ After: bending starts at the tangent modulus load, successive deflected equi- 
librium configurations must be established for which the increased column load - 
and increased internal bending are in ~The 

across the cross section of the column that changes dinae with each load in- 
_ crement, The calculation of a sequence of equilibrium positions for a succes- 
sion of small load | increments is essential, because over an appreciable por- 
= of the | column cross section the material experiences, first, an increase 
in stress under a continually changing tangent modulus, followed by a regres- 
- sion of ‘Strain that produces a stress reduction as determined by t the elastic 
The successive equilibrium configurations must be determined very ery pre-- 
= for a sequence of very small increments of unit rotation, with consid- r 
eration of the values of tangent modulus, It is not ‘prac- 


It is necessary to use a expreasion that will simulate consistently 
both the tangent modulus and compressive stress as a function of strain to as" 
great a numerical precision as may be required to give consistent numerical | 
results, | _ Duberg and Wilder? used the generalized stress-strain curves of W. 

_ Ramberg and W. R. Osgood.10 These provide a wide variety of shapes with — 
the added advantage thateach curve throughout the entire range is represented ’ 
by: a single expression. Ramberg-Osgood curves could readily be adapted 
to the procedure employed herein, but it was desired to take a very close al 

the behavior ‘near | the transition from to inelastic 


tee ‘The simulated properties for computer analysis correspond closely to - 


fer aluminum alloy 2014- T6. The minimum properties are shown by dashed 
lines in Fig. 2, and the simulated properties represented by various empirical _ 
_ relationships within th the ranges indicated are shown | on the same figure. The 

equations for stress and tangent modulus * within the r range 0,0032 <e <0. 0062 


| - 


14, 08 + +6200 +4, 34175 s sin 


do ames 


: ET “ae = 6200 + 4400 + cos 


“Numerical Procedure for Computing ‘Moments, and Buckling Loads, 
_by N. M. Newmark, Transactions, ASCE, Vol. 108, 1943, p. 1161. nine teil. a 


10 “Description of Stress-Strain Curves by Three Parameters,’ Ww. and 


WeR. Osgood, NACA T. N. No. 902, July, 1943. 
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FIG. 2.—INELASTIC STRESS, STRAIN, AND TANGENT MODU- . 
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‘It is s assumed that - for § stresses upto 33.92 ksi that the material is elastic ‘ote 
a an elastic modulus of 10,600 kips per in squared. For € > 0.0062, additional 
for and ET are indicated on Fig. 2. hers » 
_ The buckling model tobe considered herein is shown in Fig. 3. The central — 
seit rigid segment is of square cross section with breadth H and length A, 
7 q bounded on each end dbyr rigid s¢ segments of f equal length | “B. a Tt is assumed that 


(B+ 
B+. 


ie 4 for this buckling model seommatiete to t the Euler buckling stress for 

a column of uniform cross section. Eq. 4 is an approximation and, if B =0, in 
which case the buckling model becomes a column of uniform naquare cross. sec- a 


tion with the length A A, the e buckling stress is 


EE: Fr . 5 is written in terms of the radius of gyration of the cross siadiii R: -. 
~#H/v12 , the following approximation of seal Euler buckling load for a square col 


which L is as the length of the column. 

con Although the buckling load by Eq. 6is 2.7% less than the correct value by 

Eq. 7, the _ approximation in the case of the buckling model will improve, and 

the error grow less, ,as the length of the rigid segments increases, If the elas- 

tic modulus E in Eq. 4 is s replaced by the > tangent modulus ET (corresponding 

to the identical stress level), the formula will correspond to the tangent mod- 

ulus formula for column buckling. Similarly, if the reduced double modulus | 
Z 7: given by Eq. 1 is substituted in Eq. 4, the reduced modulus buckling stress will | 


be obtained. calculating En, th the value of ET corresponds, of course, to the 


| 
ie. _ which a simple equilibrium evaluation in the deflected position, using the mo- &g 
bo 4 ment area procedure, gives the following buckling stress in the elastic range: | 1 i, 
: 
— 


i If additional load above the tangent modulus load is now applied to the “ in-— 
elastic buckling model” or “ strut,” as it will be denoted herein, the strut will: 
_ ‘start to bend. This bending will cause the stress on one side, as shown, to in- 
 erease by Moy. The subscript denotes this as the first increment of stress 
above the tangent modulus value. If 4oj were thought of as an infinitesimal | 
_ quantity do, in the limit, just as bending was initiated, in the limit, there would | 
be no change in column load, and the tangent modulus would apply over the en-— 
tire c cross section. There would be no strain regression anywhere, and AC}, 2 
as shown in the figure, would be zero, If, however, Aoq is an actual finite dif- 
ference in stress level, there will be a finite increase in column load above - 
‘the tangent modulus load. if the tangent modulus governed the stress-strain 
relation over the entire cross section, any increase in load above the tangent — 
modulus load would be supercritical; the tangent modulus load itself would be 
the maximum buckling load capacity, if the tangent modulus were to govern the _ 
7 stress-strain relationship over the entire cross section. Thus, if a finite Ad _ 
_is thought of as representing the first increment of unit curvature, the strut 
cross section at the point of maximum deflection will have to develop enough Y @ 
; _ bending moment resistance tobe in equilibrium with the resultant moment due | 
_ to the external load P, increased by MP1. This increase in bending resistance 7 
can occur only if some strain reversal takes place. Thus, if the tangent mod- 4 
ulus stress is thought of as a new reference for stress changes, there will be. 
in effect a movement of the “neutral axis” inward from the edge of the cross > | 
2 section by an amount AC1,: and there will denen stress regression on this convex — 
side of the column : as indicated by 
_ If a second increment of rotation is permitted , corresponding toa further 
increase in axial column load and increase in internal resisting moment, there | 
will be a further movement during this second increment somewhat smaller 
than the first one, as indicated in Fig. 4 by ACg. This is the only way in which © 


ternal moment caused by the load P + AP} + AP2. . The increase in compres- 
sive stress indicated by Agog will be less than Ao] in the case of the typical in- 
elastic stress-strain curve in which the second derivative of o with respect to 
€ is everywhere negative, In conducting calculations of this type, it will be 
a that the average > tangent modulus during each increment of rotation 
_ governs the increase in stress that occurs during that particular | rotation and — 
sit is to be noted that within the region “C, ” the stress will first increase > 
according to the local ‘modulus and thereafter, decrease according to 
Now the change in rotation and its axial load, moment, and “Cc” for 
; — Nth increment of rotation must be considered; it is shown in Fig. 5 as the 
difference between the strain distribution at increment N and at the preceding _ 
increment N- 1, An appreciable variation in tangent modulus will be realized 
= the cross section of the strut. _ However, the increments of load and mo- 
ment, during any particular rotation increment, can be closely approximated © 
if the tangent | is determined by a single value tant decreases with each 
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= 
"stress varied ‘the increasing range, this would be the 
cation at which the total force increment could be considered as concentrated, — 
and this will give a close approximation for the meeesmaatinn distribution that ac-_ 


"4 Thus, referring to Fig. 5, the strain that will be used as an meal of of the 
changing tangent will be determined by 
€ cad 
actual ‘strain distribution entire cross section will be assumed to 
vary linearly, and this index strain, given by Eq. 8, is merely for the purpose 
- of evaluating a tangent modulus for which an assumed triangular distribution 7 
-_ of increased compressive stress is assumed to represent satisfactorily the 
actual non-linear distribution of increased compressive stress for that par-— 
ticular load. ‘The foregoing explanation | is illustrated graphically in Fig. 6, 
7 which shows a possible stress distribution at a particular bent equilibrium po- 
* sition as a solid line. In the analysis, the curved solid line portion in the right | 
portion of the stress distribution | is replaced the dashed- which inter-_ 


Referring again to Fig. 5, the magnitude of loads, moment, and *C” will i 
determined after the Nth Ao that is introducedas a result of f bending that com> 
mences at the tangent modulus load. 


On the basis of the arbitrary rotation AO, as s shown in in Fig. 5, 5, the following 


equations may be written for and AMN: 
-2 2n [( (H- 1-Cy- 1) 


acy - (E-E ETN) 


Each during a particular nume — 
— 


BUCKLING 


2 (H+2 Cy- 1) ETN + (3 H-2 Cn-1)E 
1 Op: (E- +3.N? (H-2 Cy-1) 1) ) Eqn) aC 
E-Bay) ECR] 


Int Eq. 13 values | of CN-1 will have been determined in the preceding step, 


same is true for Eq. Thus, the unknown quantity in Eqs. 13 and 14 is 
ACN which appears in first, ‘second and third p powers, For very small values 
of N, when the strut just starts to bend, ACN may be relatively large and, as 7 
bending proceeds, the quantity, NACN, increases. It will be noted that in every 
term) in which ACN appears it is multiplied by N. 
Eqs. 13 and 14 are combined into a single equation in terms of ACN by 
means of the iain n equation for the Nth equilibrium position which h may 


NAGA (B+ A) 


flection at the center line of the column, that de- 
termines ACN is also in Fig. 7. 


a Typical results of certain of the simulated tests conducted by means of the 5 
computer solution will now be presented and examined. 

In the numerical example used in the simulated tests, dimension A ( Fig. 3) - 

was held constant at 2 in. and dimension H at 1 in, Thus, the reduced section 7 

has an area of 1 | sq in., and the total strut load “Pp” is always identical with the 7 
‘eee: stress on the "reduced section. ‘Dimension B was varied in increments © 

to the value that would result in elastic buckling. 

_ Stress Distribution.—The stress distributions across the breadth of the re- 
duced section of the strut at increasing loads are shown in Figs. 8, 9, and 10, 
_ typical of short, intermediate, and long strut behavior. In Fig. 8 the space be- 
tween each of the stress distribution lines covers ten different equilibrium 


evaluations at increasing loads. ‘Thus, the used | in | the deter- 


- ae Eq. 15 the quantity in brackets to the left of the equal : sign is the maximum 
de 


Bd of 0. 0001 radians, As an indication of the fact that Ad was taken sufficiently - { 
small in the computer input, the following tabulation shows results for 


lm 
Gg 
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2 
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Calculated maximum or Shanley load 
in per square inch 


-4570 
0.000005 98. 9.4566 


Itis obvious t that range of Ad@’s that w were actually used indetermining— 


maximum loads the cited value of 0. 00001 gave er accuracy. ae 


|-00002 0.0015 3700 
2-00004 00029 3120 


3- (00044 2630 


| -0.00104 00074 2 2020 


Ue 


There is a a marked re- 

duction in the index tangent modulus This demonstrates the necessity of 

v considering the progressive change in tangent modulus as bending proceeds» , 
above the tangent modulus load. - Also to be noted is the very small deflection 


comparison with the column length at which is 


— 

§2 3 
2 
Gf" 
44 
acainto Fico & the oradnual inward movement of C ic noted ro 
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maximum lateral deflection for this strutat the maximum — is less than one * 
thousandth of the length of the strut. be Although the constriction in the chosen 7 
model tends to accentuate these effects, the model is most nearly s similar to , 
an nocaregii column in the short length range where the effects are most pro- f 


_ ‘Fig. 9 presents the simulated test stress distributions for an intermediate _ 
length of a strut. The inward extent of region C is ican than in Fig. 8. “The: 
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Toad but is than one thousandth of the strut ‘The variation 
Ef during bending above the tangent mo modulus load is is large, but not so large as 
‘Fig. 10 is fora relatively long strut that buckles just above the } proportional - 
The tangent modulus at this load is is ‘nearly the same as for the elastic 


ef J 
— 
— 4 
| 
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‘than 80% of the way into the center of the column and that the deflection at max- 
imum load is less in proportion to the length of the column than in the previous 
case. Obviously, in the elastic range there will be | noload increment, as buck- 
ling will occur at constant load, and C will not move inward gradually but will ; 
be QSatalltimes, 
Load Deflection Curves. —Fig. 11 shows typical load versus lateral deflec-_ 
tion: curves plotted from the tangent modulus load out to the maximum or Shan- 
ley load for eight simulated tests of different strut lengths. The dashed line 
7 shows the limit of lateral deflections out to the maximum load. The deflection 


is zero at the proportional limit, below which elastic buckling will occur, and 


increases to a maximum in the intermediate column range. OS 
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FIG, 10. —STRESS DISTRIBUTIONS Or<o M 46". 


- _ For a particular length L = 30 in., simulated tests ‘were made in which the 


strut was held straight above the tangent modulus to various stress levels and — 


i q _ then permitted to start bending. Results are shown in Fig. 12. The suggestion — 


_ for this type of test dates from von 1 Karman’ Ss agers of the second d Shanley | 


construct sequences of positions starting from 
load between the two limiting values ‘corresponding to the tangent 
_ modulus and the reduced moduli. The inclination of the ss 
_ representing the load as a function of the deflection is steepest for the 
line starting from the lower limiting load and becomes zero for the line 
— from the upper limiting load (i. e. the reduced modulus load). _ 
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Equilibrium lines have anevelope that starts from the lower limiting load 
and--at leastas long asthe stress strain curve canbe considered straight 
the deflection small- ~approaches a: asymptotically the | load computed 


totic approach to the reduced modulus load have noreal validity even though they - 


Or (Double Modulus Load ) 


MAXIMUM ‘LOAD INTERCEPT 


‘used inthis it is ‘obvious that the comments the asymp- 
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FIG. 12. - 12,—LATERAL DEFLECTION VERSUS AXIAL LOAD WITH STRUT HELD | 


modulus load. stress-strain curve were straight, of 
_ course, such a fictitious reduced ‘modulus load in the present case would be 
much greater than as calculated, as shown, at approximately 47.4 ksi. At 47. 4 7 
ksi the tangent modulus that determines the correct reduced modulus has sub- 
stantially decreased from the > tangent modulus at the tangent modulus load o 


- about 45.3 ksi. But the slightest deviation from straightness has a marked ef- 


fect on inelastic behavior. 
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To demonstrate the different load- deflection that results when the 
stress- strain curve remains straight above the tangent modulus (ET = const, ) 
load, — 13 should be —* with — 12, In Fig. 13 the curve of Fig. 12 


LOAD DEFLECTION CURVE 
Gone A MATERIAL WITH © 
CONSTANT Er ABOVE or 


CURVE 
(ALSO SHOWN ON FIG. 2) 


(TANGENT M MODULUS LOAD) 
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4 LATERAL DEFLECTION IN “INCHES 


: FIG. 13. _ LATERAL DEFLECTION VERSUS AXIAL LOAD FOR ASSUMPTION 01 oF 
‘that, determines the maximum (Shanley) load load (bending initiated at t op) i is re- 
drawn to a new scale. The pseudo oR based on ET at the tangent modulus load 
: is indicated and is seen to be greatly in error and much too large. Although 
- not shown herein, this error progressively increases as the strut is made 
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shorter. The illite between behavior a assuming ng constant ET and pre- 
- dicted actual behavior leading to the Shanley load is graphically demonstrated. 
Whereas, with constant ET, the lateral deflection and ‘corresponding column 
_ load would both theoretically increase as much as the geometric limitations | 
permit, the actual load differs but little from the | tangent modulus load, and the 


load is reached at extremely small lateral deflection, = 
Column Strength Curve.—Fig. 14 shows the curve of ‘strut length versus 
a average stress (or load) at the various critical loads for the buckling model. : 
_ These include the reduced modulus load, the maximum or Shanley load, and the — 
tangent modulus load. Also, for comparative the reduced — 


— 
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(Max.Lood) 
= (Tangent Mod Load 
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FIG, 14. —“COLUMN” -~ STRENGTH CURVE: 


and tangent modulus are given a 1 in. 1 in. column of constant 
—_— > curves tend to join, as they should, at small lengths, Although | 
the maximum or Shanley load is closer to the tangent modulus load than the 
- _ reduced modulus load, no particular conclusion can be drawn from this relation — 
‘because these simulated tests pertain to one particular material. | The effort 
= ae been so much to draw general conclusions, but rather to develop im- 


__LENGTH , L on 
} 
Although the simulated tests descri erein have pertained to an inelastic 
buckling model with rigid bars adjacent to the center bending portion, it would 


be possible | with ‘minor changes inthe computer program to simulate somewhat 
more nearly actual column behavior by some arbitrary assumption as to the 
distribution of M/(E I) along the column, This not done because it was” 
considered better to present a closely approximate solution of the model ex- 
amined herein rather than an incorrect solution of a more realistic represen- 


tation of a column, _ However, the present work is simply one “step along the 
toward the projected solution of columns with variable cross section in the 
inelastic range with correct evaluation of the buckled shape at the variousload 
levels above the tangent modulus load. In such a study, errors in any assumed 
_ deflection configuration would be reduced by an iterative procedure making use 
_ of Newmark’s9 numerical method as adapted to inelastic buckling, © 
ia It has been the aim of the present paper simply to improve understanding of 
‘buckling behavior in the inelastic range and to clarify : some concepts that re- 
late to Shanley’s important contributions. | On the basis of the simulated tests . 


conducted in this study and described herein tl the following may be 


made. These statements are » restricted to a material for 1 which ae is contin- : 


- ‘uous and negative above the proportional limit, 


re 


_-« 7 The inward movement of the termination point of strain gression gon 
_ the convex side of the column, immediately above the tangent modulus load, as 
_ predicted by Shanley2 in his Fig. 7, has been reproduced quantitatively fora 
"specific aluminum alloy within a localized rectangular cross sectional con- 
The procedure i: is Simple in concept, | the computer program can b be 


which suitable mathematical expressions may te written. 
3 % The procedure presented herein is projected as the initial step toward 
the accurate evaluation of maximum inelastic buckling loads for columns of 
variable cross-section and arbitrary stress-strain properties, = 
4, In determining inelastic buckling equilibrium configurations and maxi- 
-mum loads above the tangent modulus load in members made of aluminum alloy, — > 
‘it is essential to consider the continuing decrease of tangent modulus in the 


tte if Paes stability is involved in the case of a material such as struc- - 
turalaluminumalloy, 
If a column is constrained to remain straight inthe inelastic range, above 
the tangent modulus load, it will reach a maximum load that is greater than the 
= ‘Shanley load and less than the reduced modulus load, The deflection at the 
maximum anal will be progressively less as the reduced modulus load is ap- 
proached, 
7, The maximum or Shanley load for an ideal column of a typical structural 
aluminum alloy inthe inelastic range is reachedat relatively small deflections 
relative to the column breadth, 
nen 8. Further validation is given tothe significance of the tangent modulus load _ 
asa proper | basis for the Caen of | column design formulas and the double 


q 
| 
— 
= 
= 
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or reduced modulus load is seen tobe of nopractical significance for materials" 
“typified by the stress-strain curve for structural aluminum alloy used herein, 
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APPENDIX.—NOTATION 
The’ symbols a adopted for use in this paper are listed for co convenience 


= length of reduced section in buckling anniek; 
= length of rigid segment t of model; 


= distance from convex face of —_ segment to location where strain 


is neither increasing o1 or decreasing with increasing load, Nth 1equi- 
librium position; lo « 


= elastic stress-strain modulus; 


= reduced modulus; 
= 
@ ETN = effective sean modulus near the concave side of the bent “strut at 
the | Nth increment of load | during bending; 
= breadth of of reduced section in buckling m sania 


= maximum bending moment in bent segment at Nth increment of load 


NO. = subscript integer denoting a particular load increment; 


PN external column load on ‘Strut after load increment bend- 


= Nth load during bending; 
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= increment of elastic stress regression at convex face during Nth load 


increment of stress increase at concave face during Nth increment 


ofload; 


perenne increment of unit rotation; 
maximum deflection of strut centerline; 


normal strain; 


= compressive strain at tangent modulus load; _ 


normal stress; 
critical stress; : = 


compressive stress at ‘tangent 


Soy 
oR r odulus load; and 
om 
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‘By Philip G. Drazin, 1, and Louis N. Howard? 


Techniques previously used by writers’ for the of stability in a 
‘homogeneous flow are applied to problems with density stratification, Results 
are restricted to instabilities of large wavelength but are given in a fairly ex- 

plicit form. The effect of small density stratification in modifying the differ- ~ 
; ent types of instability of parallel flows of the shear-layer and jet typesis con- 
sidered, and simple spree formulas for the long-wave stability charac- 


4 INTRODUCTION 


a om The stability ¢ aa a two- dimensional horizontal parallel flow of a stratified 
inviscid fluid under gravity depends principally on three factors: The basic 
velocity field, the basic density ‘field, and the presence of boundaries, The 
density field also enters the problem in two different ways, through its effect — 


on the inertia of the fluid and through the buoyancy force. Because emphasis, — 


field is infinite and 1 no 10 boundaries 2 are present. In addition, for simplicity, the 
effect of density variation on the inertia of the fluid will be neglected. This 
effect 't is agood approximation if the total ‘change | of ae is small compared 
_ Note.—Discussion open until May 1, 1962. To extend the closing date one month, a — 
_ written request must be filed with the Executive Secretary, ASCE. This paper is part — 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the | 
American Society of Civil Engineers, Vol. 87, No. EM6, December, 1961. 
1 Dept. of Math., The Univ., Bristol, England, for merly Meteorology Dept. , Massachu-_ 
setts Inst. of Tech., Cambridge, Mass. 
8 Math. Dept., Massachusetts Inst. of Tech., Cambridge, Mass. : 
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with the mean density, frequently the most interesting — Poe The 


ethods presented here, however, can . also be used without this approximation, 


at the expense of somewhat more complicated formulas. 
a The majority of the studies of stability in . stratified flow that have been > 
made since the first extensive investigations in this field3,4 have dealt with 


specific velocity and density fields. Taylor, in 1931, , examined3 the case of a 
shear flow with constant velocity gradient in a fluid with exponential density — 
_ field, occupying a half-plane. This example was chosen in an attempt to model 

7 the stabilizing effect of cooling at the ground in 1 the evening on wind- -induced — 

turbulence in the atmosphere, . He found that there were no wave- like pertur- 

- bations, either stable or unstable, when the Richardson number lies between 
and 1/4, The Richardson number is definedas 


in which p and w are the basic density and velocity tie, and J is constant in 
the present case, When the Richardson number exceeds 1/4, only ae 
_ progressive waves are found, but he did not exclude the possibility oc 
damped or self-excited waves. However , subsequent studies of this flow from — : 
the point of view of an initial value problem5,6 7 have shown that it is stable 
for all J >0, although neutrally stable waves do not exist for 0< J< 1/4, as 
found by Taylor. 3 This flow is unstable for J <0, as is tot be e expected b because 
Goldstein, in 1931, examined* a somewhat similar flow, in which w' and 
§ p /p were constant in afinite interval, and zero outside it, w and p being con- _ 


: instability for J <1/4, in which the value of J is determined by the density and 
_ velocity gradients inthe ‘shear layer, Note that this differs from the flow stud- 
ied by Taylor, Eliassen, and Case in being actually unstable for 0<J <1/4, 
die, the fact that the density is constant outside the finite shear layer 
in Goldstein’ case, but everywhere decreasing upward in ‘Taylor’s s, associ-— 
with the greater stability of the latter case, 
Taylor ann Goldstein also considered several examples of flows with cor 
a stant shear, with layers of different constant densities, Other studies have — 
considered various flows with piecewise constant velocity and density fields. 


Several examples of this sort are available.8 


a “Effect « of Variation in Density on the Stability of Superposed Streams of Fluids. - 
we, a Taylor, Proceedings, Royal Soc. of London, Series A, Vol. 132, 1932, p. 499. 
— £40On the Stability of Superposed Streams of Fluids of Different Densities, ” by S. 
Goldstein, Proceedings, Royal Soc. of London, Series A, Vol. 132 , 1932, p. 524, earl 

3 *Two- Dimensional Perturbations ofa Flow with Constant Shear of a Stratified Fluid,” 


by Eliassen, et al., Inst. for Weather and Climate Research, ee ee Academy of 


© “Stability of an Idealized Atmosphere. I. ‘Discussion of Results,” by K. M. Case, 

Physics of Fluids, Vol. 13, March, April, 1960, p. 149. 

7 prik. Mat. Mech., L. A. Dikii, Vol. 24, No. 2, 1960, p. 240. English ndiitien 


8 “Physikalische "by Bjerknes, Bjerknes, Solberg, and Bergeron, Spr ing- 
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ta STABILITY FLUID 
a exam as give 
by J. Drazin in 1958.9 He studied the case 
constant (- « and found that with 


g2 y-2 

is is the n minimum value of the number 

the flow was stable for J > 1/4, unstable for some wavelengths for. 0 <1/4, 

he © found id neutrally stable waves of f wavelength 2 a-lq if 


: 


Miles has unified this picture somewhat by proving10 in general, 
a 


condition for stability is iy) > 1/4 everywhere. The conditien is not, , how- 
ever, necessary. This result is to be published.10,11 Miles has emphasized _ 
a ‘difficulty in 1 the interpretation of mney studies of this kind, ‘This: is that at- 


is, ‘it gives a 

is not necessarily always the « case, One need only recall the case of : a poem 
_ distribution decreasing upward, with no basic flow at all, where there are many 
neutrally stable oscillations (internal waves) and no unstable ones, to realize 


_ importance ce of this « objection and the consequent importance of the study o of 
unstable waves, i 


the present investigation the problem will be approached as follows: 
_—" as mentioned previously, the effect of boundaries is suppressed by y con-— 
sidering an infinite fluid; second, an attempt will be made to obtain results 
_ without assuming a specific form for the velocity and density fields, although | 
7 it will always be assumed that the density is statically stable, that is, it is de- 
te creasing upwards. Third, there will be particular interest in actually unstable 
: motions. In this connection the writers ” point of view is to regard results for a 


of a 1 homogeneous flow caused by ‘the presence of density variation, than as 
_ modifications of internal waves in a stratified fluidat rest caused by the pres- 
ence of a velocity field. Both points of view are of course legitimate as ways 
of relating results on inhomogeneous stability to more familiar situations, but © 
the former tends to emphasize the possibilities for instability. This point of 
view is encouraged by the the ' nature of the mathematical technique to to be used, 


9 Fluid Mechanics,” J. Drazin, No. 214, 1958. 

10 “On the Stability of Homogeneous Shear Flows,” by J. Ww. Miles, Journal of Fluid of Fluid Fluid 
IT «Note on a Paper of John Ww. Miles, ” by L. N. Howard, , Jour — of Fluid Mechanics, 
Vol.11,1961. 
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sufficient study of the relation of these to actual instability. Thus,although in 
3 _-—s«|_s many cases a neutrally stable wave motion is adjacent to an unstable one, that : 
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December, 1961 
which is a modification ¢ ‘of methods used by the writers for the study | of homo- - 
geneous stability problems, 

s ‘inally, attention will be restricted to the stability characteristics of long — 

waves. This restriction is at the. same time the key tothis mathematical tech- 

- nique and the principal limitation on its results. It appears likely, however, — 
that some such restriction is necessary if results areto be obtained for ows - 
with more or less arbitrary density and velocity profiles. _ There is good rea- 

sonto believe that some features of stability problems depend very sensitively 
on the details of the basic velocity and density profiles; this does not seem to © 

4 be the case in the long v wave region. ‘It should perhaps be stated explicitly that | 

7 all the results, like most other theoretical results in this field, ar are subject to — 
the basic limitations of linearized inviscid stability theory. In particular, the © 
possibility that some flows may be stable to infinitesimal ‘ciated un- 


_ stable to finite ones, or vice versa, , cannot be | examined, | 


>. 
vi MATHEMATICAL L FORMULATION OF THE PROBLEM 


‘The approach. in 1 stability ‘problems of examining the g1 growth properties 
of infinitesimal wave-like perturbations, by assuming a perturbation stream 
“function of theform @ (y) expia (x -ct). Here y is the vertical coordinate, 
ns positive upwards, x the horizontal coordinate, a the wave number of the 


the complex wave } epeed, _Theflow is unstable to waves of 


asic 


-c the differential equation can has writ- 


7 


as the amplitude of the displacement, due “al the perturbation, of the stream- 
Tine which, in the undisturbed state, is at ty. If the effect of f density v variation — 
on the inertia of the fluid is neglected, Eq. 5 becomes _ : 


The following will be made concerningthe velocity a and 


ve velocity field, w, approaches limits a as y -+ these limits" 


i are the same, they may be taken to be zero by changing to a coordinate system 
moving with this common n velocity. This is referred toas a “jet- -type” velocity 
7 profile. If the limits are ‘different it can be assumed that they are negatives © 
a each other, by changing to a coordinate system moving with average velocity | 

a the limits. This is denoted as a profile of “shear-layer” type. Typical jet 


shear layer velocity profiles are shown in Fig. 
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Yl 9 _ zontal with velocity w(y) in the positive x direction and the basic density field eke 
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STABILITY IN FLUID | 
2, It is assumed that t p is monotonically de decreasing upwards and: approaches 7 
“finite limits as y~-te« & To be consistent with the use of Eq. 6 instead of Eq. 
we also suppose that p(-@) - p(+«) is + oo), 


* typical density profile is shown in Fig, 2. 


SHEAR 
LAYER 


- L _ Vv will be typical length and speed scales of the velocity profile and 
will be sneer to define dimensionless length, TC and time variables in the 


IG. 1 1. PICAL VELOCITY FIG. 
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Assuming a relatively small ‘total change ¢ of, Pp, Eqs. 8a and 9a are approxi- 


) = (+0) 2 ply) 


giving a somewhat more e familiar form of the Richardson es J. It should 
be n noted that A(+o) = + 1 . The relation of this over-all Richardson number 


and the local J(y) = -2 is, in terms of the ‘present notation, 


The mathematical problem is. the determination of the ‘conditions under 

. ieee Eq.7 has a non-trivial solution that approaches zero att, Suchaso- 
lution will usually be possible only if c anda are properly related, The ap- 
proach to finding this “eigenvalue relation” between c anda is first to find, 
essentially by means of an expansion in powers of a for fixed c, a solution oy 

va that satisfies the boundary condition ¢1(- ©) = 0. Similarly, ¢2 is found with 

7 9(+ eo) i. If now c is to have the correct value for each a to give an eigen- 

function: Ky, = Ke for some constants Kj and Kg, so that 

“Ke $5(0) . ° 


64° (0) = Kg 


Because @ is not zero, Ky and Kg cannot ‘zero an therefore, 


Oya and O25 , of course, depend on and c, and this condition then yields” the 
eigenvalue relation, The problem is, thus, Teduced to finding some representa- = 
for the functions $1 and $2. A As -- the differential Eq. 7 becomes _ 

because it has been assumed that w ~ constant and = 0, Thus, 

_ the asymptotic form of the solution ¢1 which goesto zero at = can be taken 
to be py~ e®Y, This shows that an attempt cannot be made to find $4 by 2 
BS series in a, because the convergence of the series must become in- 
worse -- ©, However, it turns out that such a series 


SKETCH OF METHOD OF SOLUTION 
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_ STABILITY IN FLUID 
and then expands 6. Similarly, ¢9 is found setting 
2 


and expanding x in powers of a. In this | way a . representation of 4 is obtain ned 
which is valid except near y=t+o anda | representation of $9 valid except near 


ye-+e . It is because of the form of these representations that the eigenvalue 
relation is obtained by using thetwo functions and matching them in t the middle | 
s for example, y = 


0, rather than, for instance, trying to make ¢1 zero at + 
The representation of fails to be usable at+ 


In the present problem, we must also deal with the term involving J. As 
‘mentioned in the introduction, our point of view in this paper is to look at the 

problem as a modification due to density variation of the homogeneous stability 
"problem and, consequently, we regard J as small and actually find 6 and x in 


the form of double series in power of a al and J, settling, for instance, 


. 20 + @ 5 + J* +. 
Introducing @ and x instead of @ in Eq. 7, expanding in powers ers of a and J, 
and using the normalizations @/W as 


- cand x/W ~lasy the 


— ’ 
ail 
| 
(w - + f¥ ay’ f¥ (1 te) 


ry (1+ (1+ 
| te. 


(17h) 


If these equations are now used in Eq. 12, the neice up to the second order © 
in a and J are written down explicitly, the eigenvalue relation |, Lae 


- 7 Notice that this is ina relatively explicit form in that it involves only inte- 
 grals of functions which are completely known (in terms of the functions A and 
w) except for the constantc, 


COMPARISON WITH FAMILIAR RESULTS 


a simplest example of a non- homogeneous stability problem is probably a 
that the ‘Stability of the interface /betweentwo paralleluniform streams of 


ig 
-a)| [a (w2 - + (1 +a)| 2 
— 


STABILITY IN FLUID 
different densities, The solution of this problem dates from th the work of Helm- ; 


holtz12 and the in the dimensional form is: 


oo) + + oo) 


or, in terms of the a ae variables introduced previously, a and assum- ; 
ing that the change of density small, 


_ This case does not strictly fit intothe scheme used because it has been as- _ 
: sumed that the density field is differentiable. However, because the final for- 
‘mula, Eq. 18, does not involve derivatives of A it is of interest to try it on this 
problem, In fact, Eq. 18 succeeds it 


so that the second order terms vanish, and 


we =| (1 - + = a(1 + (22) 


"Thus, Eq. 18 gives” Ea. 20 exactly, at least to the order of ‘approximation in- 
_ dicated. This is of course not - surprising, for the Helmholtz flow, any wave- 
It will be noticed that w-c appears in the denominator of the integral in 
: Eq. 18, which, thus, becomes difficult to interpret if c is real and in the range 
of w.) For this reason Eq. 18 will, for the most part, be restricted to the un- 
_ stable case, in which c is not real. : There is one case, however, in which this — 
difficulty does not occur, namely in which w = 0, and Eq. 18 then gives a for-_ 


for — of waves ina stratified fluid at rest, namely: 


The first approximation J/a is easily verified to be the for- 
mula for internal waves on a —_ discontinuity, which thus holds in any cone case 
ives 


for the internal wave with c2 ~ as as a-0. Of course witha ‘continuously 

_ stratified fluid there are an infinite number of different internal wave modes. 
_ The one just described is the fastest of these. To describe the others, more 

terms a are needed, although Eq. 23 will also give the second as it stands. 


” by H. Lan mb, Cambr idge Univ. P “e s, 6t 
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December, 1961 
This can be c checked, for instance, by examining the case ase of a three layer sys-_ 


STABILITY BOUNDARIES IN THE SHEAR-LAYER CASE 
A velocity profile of the shear-layer type should no now be considered with 
@gsmalland J small, even fairly small a. . The Eq. 18 has: the 
— 
As a result the flow is unstable, the instability being essentially a weakened 
he form of the instability of a shear layer in the absence of density gradients. So 
long as J/a is smaller ‘than. c is, in 1 first. ‘approximation, imaginary, 
is no difficulty about convergence of the integral. By using the first ap- 
_ proximation to c in the integral, relatively simple formulas for the second ap-_ 
_ proximation can be obtained. However, it is natural to look for the boundary of a 


this regionof instability: howlarge can J be made, for agiven a, beforethe in- 
stability is suppressed? The first approximation of course gives for the sta- a 
bility boundary the equation = However, , because this value of J makes 


a ond approximation are not defined if this value of c is used in them; conse- 7 

quently, caution is required in investigating the stability boundary. _ Clearly, — 
what is needed is the limiting value of the integral as c becomes real, that is, — 
asthe stability boundary is ageronched. This can be obtained, but the analytical 


details are somewhat lengthy. The result will be given in the simplest | case SW 
in which w is an odd function with > 0 and A is odd, In this case, ap- 


— zero through pure imaginary © yaioee, and the second approximation 


7 


ole 


A Eq. 27 gives a remarkably simple representation of the stability , boundary. 

_ for essentially arbitrary odd profiles of the shear layer type. Examples will | 

~Serve to illustrate its usefulness, ‘The is ‘essentially a a problem solved 
by in 1931: 


| — 
: 
4 
3 
Solving this for 
— 
— ~ 


ot 


r ca th 


_— Eq. 27 gives a fairly good approximation for a smaller than about 
1/3. If the quadratic approximation is taken literally, it is found that J has a — 
maximum of 3/8 at a = 3/4. Now for this problem the local Richardson num- 

‘ber J(y ) in (-1,1) and is not defined for |y| > > had condi- 


tion for ‘stability (12 > })e can be replaced by the c condition “ - p an 


in this form it can seen his” 


condition | gives stability. for a = q in this case, Thus, even though a. 


‘much too for validity of quadratic approximation, the value ob- 
tained from it does not hopelessly miss the absolute upper bound > given by 
_ _ The second example is w= As = tanh y. This problem was s solved by. J. 
Holmboe, in 1960,13 with the exact result J = a (1 - o). This appears to bea 


favorable case for of Eq. 2 27 | exact result i is 
ratic in a. ItyieldS 


=f (tanh? Le 


00 
so o tha t Eq. 27 is actually exact inthis case. | : 
is perhaps of interest to give the details of this somewhat 


lly. The local Richardson number is: SS 
ay 


at(y) 


Thus, the condition th that J J (y) b be everywhere 2 23 
“The eigenfunction and eigenvalue relation a are in nthe case 
; a corresponding to the stability boundary in which c = 0, In this case, as men- 
tioned previously, the eigenvalue relation reduces" tod j= -a), and the 
corresponding eigenfunction is 
sech® -y |tanh 
as is easily checked. may be ‘compared with that obtained 
from the formulas for the expansion coefficients 


EM STABILITY IN FLUID 
In 
Th 
result 
— 
— 
— 
J(y) = J = J cosh2 * ( 31) 
— 
4 
ae __ 13 Unpublished lecture notes by J. Holmbee, Univ. of California, Los Angeles, Calif. § 7m 
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give the results in the ¢ case of the Putting c and evaluating the the inte- 


gives, (setting T = tanh) for short) 


‘T)+2T Sin 

1-T 


a2 on 1 the stability boundary, 


Using t the: values given in 33, to order a2, yields 


a 2 
~*~ 


a = 20 e-ay T 
= a ; shows. that @ ~ 2% e-@y for y ~~, so that to make the normalization 
of agree with that used asa basis for our formulas we must re- 


‘place o by 2-2 9, Because 


ta 
us 
x... = 2(T-1) - (2741) 
&g 
im 
— 
q 


it is seen that the expe expansion presented here | agrees with the exact solution on 
the stability boundary. Because the exact solution is not knownfor c # 0, fur-_ 
ther comparisons cannot be made, However, the use of these formulas will 
be illustrated by obtaining the eigenvalue relation for Holmboe’s example in 
the unstable range J <a(1 - a), to second order in a, For this Eq. 18 will be 
“used. After introducing T = =tanhya as) of Eq. (18 becomes: 


be shown that cy = 0. This is physically fairly clear from the nainnanehes 
of the w and A fields which shows that there is no preferred directionof propa- 
gation of unstable waves, and that one would, consequently, expect the instability 
to be simply an exponential growth of the initial sinuous perturbation, — 


7 


a J j2 


on —EM6 STABILITY IN FLUID 113 
&§ 
q 
&g 
aim 
Writing Cy + and Keeping in min at It 1S assume that Cy >0,the 
ald + e%) - 25 + 4 + | 
+ 2 ic | + tan-1 |———_} |" 
“ 
J=a (1 - 2c, tan-l(— | -2 


7 ro first order, one obtains, to second order | 


From Eq. 43 it is a simple matter to compute curves of constant Cj in the 


a plane, 
‘stant growth rate; acy = ‘constant. “Such” curves shown in 1 Fig. 
the maxima of these curves one can determine the fastest growing wave num- 
_ ber: for any value value of J <1/4, ‘This is th the wave number one would 1 expect to ob a 7 


_ serve in an in experiment, — ‘The | locus of these maxima is indicated by t the curve 
will be noticed that the ge. general formula, 27 does not agree with the 


Boca of Drazin9 mentioned in the introduction, which has J or order a2 rather 
nential in the whole plane and sothe assumption that p approaches finite limits" 
JET TYPE VELOCITY 


es a single maximum n of, say, and ‘approaches 2 zero ) at - + have just two 
_ types of long wave instabilities. The nature of thesetwo modes of disturbance 
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is indicated by Rayleigh’s names for ene sinuousities a and varicos: ities. . E}. 
18 for the jet case and with J = 0 can be written Te ee 


-a 


‘The sinuous disturbance has c — 0, through nearly aad imaginary values 
a — 0, and to first 


= layer case to” stay in the unstable region. In fact, Eq. 18 shows that 7 
must be small of the order of a2 to do this, and as a first approximation to the 


equation of the stability 


‘This refers to the sinuous instability. The varicose se instability isa little 
‘more complicated. Here in the case J = 0, c +1 (the maximum of w(y)). 
_ through ¢ complex values, as a@ — 0. The precise behavior depends on the nature 
of w(y) in the neighborhood of its maximum. Inthe presence of density varia-_ 


- a eyo ‘situation prevails. Eq. 7 shows that in the limit of a -0 and 
a 0 w2 (9/w)" can be expected to be constant, and because it is zero at ec, 
ally 2 zero. This implies that 9 /W is constant throughout intervals in which — 


W # 0, but this constant value may change at a point at which W = 0. I poor 
a change occurs, 9 will still be continuous because - 0 at the e jump in the in the 


constant, but will not unless also there, 
Thus if a discontinuity in vorticity is to be avoided, c must is the 
- maximum value taken on by w, as @and J approach zero, (Only the simplest : 

0 


 jet- type case, in which w' vanishes at just one point, say y = is considered — 


here, ) Ifc approaches the maximum of w, for example, 1, however, the inte-— 
gral in Eq. 18 diverges sothat the secondorderterms are not really of second _ 


It can be shown, however, , that the higher order terms that have been» 7 


order. 
If is 


"omitted in . Eq. 18 are still of smaller order than those retained. 
equal to 1 + y, in which it is assumed that y is a small complex number with | 
argy <7 - ¢€ for some > 0, the integral in Eq. 18 can be estimated, and 


‘retaining only the largest term for small y Eq. 18 yields 


ae _ Thus the jet, although unstable to long waves, is less so than the shear lay- ; x 
inthe sensethat the corresponding growth rate is smaller. Similarly,when 
inwhich, if y = |7| 7-3/2 means |yj-3/2 
ot 


_-‘Thus if J is small enough in comparison to a, , the right- hand side of Eq. 48 
is negative; o = 27/3, and instability is obtained, with 


7 -(0))] [a- +a(o) 12/3 

_ It is clear also that the first approximation to the ‘stability boundary : is 


_ The research of the senior author was partially supported by the Office 7 
Naval Research under contract N-onr 1841(18), 
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- RANGE OF YIELD ) CONDITION IN IDEAL PLASTICITY — 


By R. R. M. 1M ASCE 


SYNOPSIS 


The properties o of ideal rigid-plastic solids are used to define bounding yield 
conditions that remain valid when test results are available for only one stress 
state, such as simple tension or simple shear, For material in which yield is 
_ insensitive to changes in mean | stress, the 1e bounding yield conditions are > shown - 
to be the maximum shear-stress criterion of Tresca and a new criterion n des 
_ ignated asthe maximum reduced stress criterion. Use of bothcriteria enables _ 
the best possible bounds tobe obtained for the carrying capacity of a structure - 
in the absence of further information, Parallel results are obtained for ‘ma- 
terial in which yield strength is alinear function of the mean stress, and some 
examples are given, 


- Note.—Discussion open until May 1, 1962. To extend the | closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. This paper is part — 
of the copyrighted Journal of the ta en Mechanics Division, — of the 


aie Plastic Setematinne, ” by A. A. Gvozdev, Proceedings, Conf. on Plastic Defor- 
_ mations, 1966. In Russian, edited by B. C. Galerkin, AN, SSSR, 1938, p. 19. on 
translation: Internatl. Journal of Mechanical Science, Vol. 1, 1960, p. 322, 


— is 
&g 
3 
| 
aes | _ The powerful variational principles of the theory of ideally plastic solids” A 
ie . have been used primarily inthe development of approximate solutions to well- : 
“4 * |. defined problems in which the yield properties of the materialare supposed to 
_ be known exactly. The carrying capacity may then be bounded from above and : 
by means of the theorems of limit analysis, developed by A.A. Gvozdev2 
i 
| 
Prof. of Enere. Science. The Univ. of Michigan. Ann Arbor. Mich. 
&g 
— 


and by D. C. Drucker, Ww. Prager,and H.J. or by the use of escribed 
i: inscribed yield c conditions. 4 On the other hand, the question of finding ap- 


proximate solutions when information about the yield condition is incomplete — 
does not seem to have been investigated, despite its technological importance. a 
a this paper, the problem of finding upper and lower bounds to the yield | 

: : ell of a body is studied for the case when the yield strength of the material is 
known exactly for only one or two stress states. An isotropic material in which 

7 [ yield is independent of the mean stress is considered first, and then the anal- 
: ysis is extended to material in which the yield criterion is linearly dependent 
onthe mean stress, The former model is widely used inthe analysis of metals, 
‘whereas interest inthe ane ‘springs from mattempts to apply | the ee of plas- 
Notation.—The letter reniaie adopted for use in this pa) paper r are defined where 
fey first appear, in the illustrations or in the text, and are arranged lalphabet- 


MATERIAL INSENSITIVE T¢ TON MEAN STRESS" 
‘For this paper, an ideally plastic, isotropic material, in which the y yieldcr cri- 
_ terion is independent of the mean stress, will be considered. In view of the 
Saati of — the orientation in space of the principal stress o_o 


principal stress components 01, 09, and 03. It is then possible to en Rr 
the yield criterion by use of a surface in the principal- -stress space having - 
‘components 01, 02, and 03 as co- -ordinates, In the case of yield criteria inde-_ 7 
pendent of the mean stre e will be by generators at 


that is, , parallel to the axis is making equal angles with the co- 0-ordinate or: 

and any int intersection of the surface with a plane represented d by Eq. 1 will be 

ee ‘Two cases of intersections of this type are shown in Fig. “3 One c can first 
consider the case in which the results of tensile and compressive tests are 
available (Fig. 1(a)). In the case of metals, tensile and compressive yield 
strengths are often nearly equal, suggesting the additional assumption that will 
+: be made herein, that reversal of the sign of a stress does not alter the stress | 
magnitude at yield. Test points are represented by the six small circles in 
Fig. 1(a), because any one of the three principal stresses can be taken as the | 
- non-zero component, Alternatively, considering a case in which the result of 
q shear test (such as simple torsion) is available, the six points indicated by 


the small circles in Fig. .1(b) will then have been determined experimentally. .- 


3 “Extended Limit Design Theorems for Continuous Media,” by D.C. Drucker, Ww. - 
Pra er and H. J. Greenberg, Guarter’y of Applied Mathematics, Vol. 9, 1952, p. 381. 
“A Theoretical Investigation of the ect of Specimen Size in the Measurement of 
” by R. Hill, Philosophical Magazine, Vol. 41,1951, p.745. 
5 “Soil Mechanics and Work-hardening 7 Theories of Plasticity, ” by D.C. - Drucker, R. 
E. Gibson, and D. J. Henkel, Proceedings, ASCE, Vol. 81,1955, p. 791. = | 
6 “Mechanics of the Triaxial Test for Soils, *by R. M. Haythornthwaite, 
ASCE, Vol. 86, No. SM 5, 1960, 
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YIELD 


‘When finding the smallest surfaces, it is sufficient to note the requirement of - 
convexity for any yield surface associated with a stable plastic material, 7 in 
_ order to arrive at the two hexagons shown. The largest surfaces are found by 


noting, in n addition, that the cross sections shown must have 30° symmetry, so 
_ that reversal of all three stress components does not influence yield. — 


eg smallest yield surfaces that can be em through the test data inthe two cases? 


@ (a) TENSILE- TEST RESU | SIMPLE SHEAR- TEST 
AVAILABLE RESULTS AVAILABLE 

- FIG, 1.- —BOUNDING YIELD CRITERIA FOR MATERIAL INSENSITIVE TO MEAN 


The hexagons marked 1 in Fig. will be recognized as representing 


familiar maximum shear stress criterion of Tresca: 


(jor - o2|, |o2 - 03|, |03 - o|) = 09 = 
in which | o9 is the yield stress in 1 simple tension, and TO is the yield stress in 
- simple shear. The hexagons marked 2 in Fig. 1 represent criteria in which a 
restriction is placed on the value | of the maximum reduced stress (which | in 
Fig. 1 is measured to the stress axes): 


-2, FO (3) 


“max - - o| 


a a) More Fundamental Approach to Plastic ‘Stress- -Strain Relations,” by D. C. 
Drucker, Proceedings, , ASME, Ast ‘S. of f Applied Mechanics, Chicago, 


1961, p.487. 
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The shown in Fig. 1 and represented analytically by Eqs. 2 and 
owe their significance to the property that all possible yield conditions must _ 
lie between them. A lower bound obtained by the use of the inscribed hexagon 

be a lower bound for all stable materials irrespective of the details of 


yield condition, and an upper bound obtained by the use of the escribed a 
be similarly an upper bound in all 


a 


> 
reentt hexagons, but did not make “this property dependent on the ae 
of tensile test data, Neither writer anticipates the particular use of the cri- 
torts that is suggested here. 
the computations are greatly simplified, if 
a ume expression can 1 be found ‘for the rate of energy dissipation in terms 
of the strain rates, The theory of the generalized plastic potential10 indicates — 
that the associated strain- “rate vectors coincide with th the outward drawn normal — 
— 


in which Xp are positive constants : and the fp are vege ular functions, the « equa- 7 


tions of various parts of the yield surface — 


cht 


in the case of the maximum shear-stress criterion, and - 
: inv which y is ‘is s the e shear : strain, in the case of the maximum reduced stress cri- 
«Still further restrictions on the range of the yield criterion can be made if 
the results of more than one type of test are available. If, for example, results” 
for both tension and torsion tests were known and these results were consistent | 
the maximum shear-stress” criterion, an examination of Fig. 1 reveals 
that the criterion must hold exactly for all other stress states. On the other oe 
hand, if these tests were consistent with the maximum reduced s stress criterion, — 7 


ar 8 “On. the Inhomogeneous Deformation of a Plastic Lamina ina Compression Test,” 
by R. ‘Hill, Philosophical Magazine, Vol. 41, 1950, p. (733.0 
"ee 9 “On the Development of a Theory of Ideal Plasticity,” by D. D. Ivlev, ne 

(AN, SSSR), Vol. 22, 1958, p. 850, (in Russian). English trans— 


lation, Pergamon Press, Inc. p.1221. 
_ 10 “Stress-strain Relations, Uniqueness and Variational Theorems for. Elastic-plas-_ 
tie Materials with a Singular Yield Surface,” by W. T. Koiter, Quarterly of Applied 7 


— 
— 
| 
— 
a3 It is ten a simple matter to show that the rate of energy dissipation is 7 he 
‘ 
i 
| - 
} 
| 
| 
— 
— 
a 
| 
— 


that criterion n would be { the o1 only nti. one. = both ca cases, the range » of pos- 

_ sible criteria is reduced to a band of zero width re 
— In other cases, it is evident that additional results from only a very few 

_ combined stress tests would be sufficient to confine the yield surface fora 

stable material to an extremely narrow range 
As might be expected, application of the yield criteria developed previously 
leads to bounds on the carrying capacity which, in some cases, are much closer | 
than can be obtained by a simple factoring process based on inscribed and 7 
_escribed yield surfaces of the same shape. If tensile-test values are available, 
an upper bound based on the Tresca yieldcriterion must use the outer hexagon 
shown in Fig. 2, whereas an upper bound based on the maximum reduced stress 
_ eriterion can use the inner hexagon, — The latter will always be at least as fa- 
vorable because, as a consequence of the lower-bound theorem of limit anal- 
ysis,3 increasing the yield strength of a body in any zone can never reduce the 

_ A Plane Stress Example.— Acircular plate of radius R and thicknesst << RB, 

d simply supported around the edge and subjected toa lateral pressure p uni- 7 
_ formly distributed over the top surface will now be considered. _ Supposing t the 
“yield stress of the material in tension tobe, given, it is desired to find the clos- - 
est upper and lower bounds to the carrying capacity. For the lower ‘bound, a : 

solution using the Tresca yield criterion will be appropriate, ‘This 
given by H. G, Hopkins a and Prager. 11 ‘The corresponding solution for re- 

"duced stress criterion will be developed herein, 

ih the analysis, the usual assumptions of plate theory will be adopted. 
the effects of shear are “neglected, and if the lateral pressure is small com- ¢@ 
pared with the stresses in the plane of the plate, then the state of the plate —_— rf 
be approximately one of plane stress, the non- zero principal stresses lying in 

the plane of the plate. _ Adopting the maximum reduced-stress criterion, the 
yield curve drawn in terms of the radial moment per unit length, Mr, and the 
circumferential moment per ‘unit length, Mg, will be as shown in Fig. 3. Be- . 7 
cause Mr is zero at the outer edge: and equal to Mg at the center (by symme- 

the state are expected tolie in the range ABC in Fig. ~The equa- 


and that of BC is 


The Load Carrying Capacities of Circular Plates,” by H. G. dW. 
Prager, Journal of the Mechanics and Physics of Solids, Vol. 2, 1953, 


| 
— 
4 
— 
; 
| 
= 
Mg -gMr=Mp (8) 
me _Eliminating the shear force from the radial-equilibrium equation by use of eS 
| 


‘FIG, 2.—ALTE 


FIG, 3.—REDUCED STRESS CRITERION F ING 


— ITIVE ESCRIBED YIELD C 
“tw 


YIELD CONDITION 


induce sag of the plate in the same direction as the action of the pressure, Pp. 
‘Substituting Eqs. 9 and 10 into 11 11 and integrating | the ) resultant first-order 


Gaim 


Substituting the conditions ‘My = 0, Mo = Mo R ont My =Mg 
= Mo at r = 0, and equating the resultant expressions for for Mr at r= "a ie R, the 


transcendental equation for is obtained 


The root in the range 12>B2 0 is B= 0.€ 6242 which, on substitution into Eqs, 


ps 6.852 M0 


= pressure from Eq. 15° is statically 2 admissible and § so is a lower bound. It 

an be shown to be the actual collapse load by associating a velocity field, 

7 ”" The velocities” are determined to within an arbitrary constant through the 

position of the state point in Fig. 3. According to the theory of plasticity for 
stresses, 12 the generalized strains corresponding to Mr, Mg will 4 

a form a vector coinciding with the outwards drawn normal tothe yield criterion, 

corresponding axes are ~The generalized strain rates cor- 


12 “General Theory of Limit W. Proceedings, 8th Internatl. Con- 


— . 
— 
—- 
— 
= 
on AB an 
— 
— 
¥; 
s 


and in a similar fashion the 1 normal to Eq. 10 leads i 


1.240 p2) w, wo 
in which Wo is the deflection rate at the center, A velocity seine has 
been associated with the he statically admissible load, Eq. 15, and the solution is" 
——s- The estimates of collapse pressure obtained, using the two criteria, can 


now be compared. In the case of Tresca’s criterion, ‘the pres- 
sure given by Hopkins and Prager!1 is oe 


| 
and this will be ‘the best lower bou 


a) in which w is the transverse velocity. The normal to the line described by Eq. Bi 
: 
(20) 
Setting the deflection as zero at the supports, and noting be 
nd. Comparing Eqs. 15 and 23 
— 


YIELD CONDITION 


and the pressure is bound to within 14. 2%. 
Inc contrast, if the Tresca yield criterion were - to. obtain the upper b bound, 


_and the pressure is bounded to within 33.3%. 
"MATERIAL SE SENSITIVE VE TO ‘MEAN STRESS 
An ideally plastic, isotropic material in which the s ei level at yield isa 
‘me function of the mean stress will nc now be considered, This constitutes a 


TENSILE- TEST RESULTS SIMPLE. SHEAR-TEST 
AVAILABLE RESULTS AVAILABLE 


— 
— YIELD FOR MATERIAL SENSITIVE E TO MEAN ST STRESS 


compressive ‘tests by the. are available; ‘Fig. 
_ the case in whichshear tests are available. The limiting lines are obtais obtained ed by 
arguments parallel to those u used previously. 

= When expressing the bounding criteria mathematically, it is convenient a 
Bessy the cohesion, c,and the angle of internal friction, ¢, in the fashion in 
which they are commonly used in statements of Coulomb’s yield criterion. Any 


line that cuts the octahedral axis can then be defined interms of c values 


a 
lee 
| 
| 
| . Ny so 
a “first step in the generalization of the foregoing analysis. There is evidence _ 2 
that it may be appropriate for certain soils and possibly for some concretes 
___Intersections of the yield surfaces with the planes determined by Eg. 1 are 
“ 
‘he 


6 


‘TABLE _EXPRESSIONS FO FOR THE YIELD 


where 2-sin 


= (s2+03)(1- -sing) /2+ccos> 


mine (01-0, G2-0, = 4 ccos$)/(3- sing) 


(1+sing) +92) /2-(1- sing)d3 = 
205 = [ (2a- ind) 03 /(142)) 


03; wes sind/sing' ) 


(0, > > 


una = — 

= 
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| 
q 
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= 
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CRITERIA ILLUSTRATED IN FIG. 4 


(1-sing' =O 


STRAIN-RATE BQUATIONS 


l- -sind) = 


YIELD 


= 02 > > 03 


=) 


(1- 


ic 
sing)ert(1+sing) (estes) = 0 
= 

Ed: cis -2(1- -sino 


)(1+sing) 
-sing! 
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02 > 93) 
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sing) €3 = 
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= 2 l+sing: sing:sing-1 


-2sing)€i+ (1+sing) (« 


> 
> 52 = 


= 4 bs 
— 
— 
tm 
sing') 
— 
ta 


= 


to the criterion that contains the line, 
‘The algebra is straightforward and will not be given herein; the resulting ex- 
pressions for the various bounding criteria are summarized in Table 1 . When 
tensile and compressive tests are available, these are not necessarily | con-— 
sistent with the same generalized Coulomb yield surface.6 Therefore, values | 
cand ¢, corresponding to the compressive test results, and c' and ¢' corre-- 
sponding to the tensile test results, are introduced. The ratio a = $/@' is lim- 
requirement of convexity to 


7 


the extreme values $ representing the attainment of triangular cross sections to 
the yield surfaces A and When a= = 1, surface A reduces to the Coulomb 
yield surface (Table 1), 
_ The flow equations in the Col. 3 of Table 1 are deetent by substituting the 
analytic expressions for the yield onetes into Eq. 5 given in Col. .. 


At point 


which also applies to al all other points ¢ on the surfaces; the apex can be + consid- a 


- ered as belonging to every side, and on any one side the p projection n of stresses 


4 in the direction normal to the side is always the same. In the case of surface 


(29), 


in which ty is the value of largest shear strain rate. Equally simple 
; ee do not appear to exist for the other surfaces, and Eq. 28 has to be 


— 
— 


iow Strain.—In the case of plane strain, all the bounding criteria intro- 
duced previously can be expressed in terms of equivalent Coulomb criteria, 
with a suitable choice of the constants. A solution that has used the Coulomb 
criterion can then be readily adapted to paore to any of the other criteria by | a 
In plane strain, one principal Strain rate is to be be set equal to zero, The 
>€9>€3 is already establishedas a consequence of the ordering» 
_ 01 2 02 2 03 adopted in Table 1, and it remains to determine which strain rate 
is to be zero, The dilatation is everywhere positive, so €14 + 0. wag* = 0, then 
7 4 and €9 are both greater than zero, but by the flow equations in Table 1 they 
are both of opposite sign in every case, which is a contradiction; hence, €9 = 0. 


For surface A, Fig. 4, substitution of €2 = 0 in the flow equations, Table 21, 
corresponding to th the two types of corner, 


+ (1+ in 9) =0 (30 


| 


4 
— 
— 
— 
— 
— 
— 
aig 


4 YIELD CONDITION — 


when oj = 02 


sin 01 - - sin ¢) 03 = = 2c cos 


( sin ¢') oy - (1 ) 03 = 2c’ cos ol 
"when o< ¢' . Thus, for plane strain, surface A reduces to Coulomb’ Ss criterion, — 


_ the constants c and ¢ being used when ¢ 2 ¢' and c' and $' when $< . 
‘For surface B, Fig. is one type corner. "Substituting the 


corresponding stress 09 = = 05 ip the yield ‘criterion, ‘Table’ 1,1 we obtain, after 
some reduction, an equivalent Coulomb yield criterion in which the effective 


values 6* and c* are defined by 


sin 


to th the Coulomb criterion when the value oe 
fe which , holds ata corner, is ‘substituted. a is as expected, because this value 
of a2 represents the the state of simples which the experimental data tawere 
For surface D, Fig. 4, substitution tee = 0 in the flow equations, Table 1, 
: corresponding to the two types of corner, ae 


= = > 03 and 


2 sin 6) € 


a > 03 and 
(1 - sin + (1+sin ¢') €3=0 (30d) 
when 01 > 02 = 03. Examination of the plane stress cross section (not shown) 
im =—=sreveails that Eq. 30(a) defines a vector lying between normals to adjacent flats 
* providing ¢ 2 ¢' and Eq. 30(b) defines a similar vector providing ¢ < ¢' Sub- 
a -_ gtituting the stress conditions associated with Eqs. 30 in the equation for the ig 
ou 
in (3 - sin 6) (1 - sin 9") 
| 
: 
— 
a+ (1+sin@) €g =O (3d) 
4 
— 
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when 01 > 02 = 03. Examination of the plane-stress cross section of the yield 
surface “shows that Eq. 34a is inadmissible because it is not the normal toa 
_ supporting plane of the yieldsurface. Substituting o2 = o3 inthe e yield criterion, 

Table 1, an equivalent Coulomb yield criterion i is obtained in which t the effective an 


values * and c** are defined by 


1-2sin¢ 
1- sin 


As an s an alternative approach, the same effective values of ¢ and a6 for substitu- 
“tion j in an equivalent Coulomb yield criterion can be found by comparing the 
corresponding flow equations given in Table 
AML problems of plane strain have now been shown to reduce to equivalent 

_ problems formulated in terms of the Coulomb yield criterion, In the case of | 

materials insensitive to mean stress, the parallel result is well known. —_— 
- ever, it must be recalled that, in the presence of dilatation, the material is no — 


ze in a state of simple shear and the argument used in the latter case land 


Plane-Strain '—Consider the identation of a half space by a long, 


- stants in Coulomb’s law, so that a = 1. As stated ' ‘by L. Prandtl,13 the inden- - 


Frigid punch i in 1 the particular case when | the tensile and compressive strengths j 


tation pressure is 


Prandtl’ s solution is complete in the sense than an incipient velocity field can 
be associated with the stresses in the deformable zone, and the material can 

be shown to be necessarily rigid in the remainder of the half space,14 15 oe 
me: In plane strain, surface A, Fig. 4, reduces tothe Coulomb criterion with teed 
constants c and while surface B requires the new constants defined 


and because the ‘resultant stresses velocities < are e consistent 


_ the various yield surfaces of Fig. 4, in the plane-strain case . The a 
2 r i and lower be bounds sare shown inF 1 Fi 5 for various ¢, the ‘Tange be between 


ichen Gesellschaft der ZU Gottingen, Mathematisch- Physikalische Kla- 


_ 14 * Mixed Boundary Value in Soil Mechanics, "by R R. T. Shield, Quarterly 


of Applied Mathematics, Vol. 11, 1953, p. 
“hp iastic Potential Theory and Prandtl Bearing Capacity Solution, ” by R. T. Shield, 
ASME, Journal of Applied Mechanics, Vol. 21, 1954, p. 193. 


: 13 “Uber die Harte plastischer Korper,” by L. Prandtl, ‘Nachrichten von der Konig- 
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5.— PUNCH INDENTATION (PLANE STRESS) FOR FOR A 
_MATERIAL SENSITIVE TO MEAN STRESS 


the ond lower curves representing the maximum » variation in 
that can occur with convex, conical yield | surface 
CONCLUSIONS 


Certain yield criteria have been established as intrinsically suitable for 


may b be, without | any reliance on a detailed | knowledge of the actual yield cri- a 
terion. All the criteria are piecewise linear, , with the attendant simplifications: 
the analytical work that are well known, 
«It is perhaps worth emphasizing that results from af few aggre may | serve to 
_ define the yield criterion exactly, 


_ of the bounding criteria. If, for example, a metal develops yield strengths ce 
tension, and torsion consistent with ‘Tresca’s 3 criterion, no further testing is 


i applies for the Coulomb criterion in the case eas ar material for — 
_ yield is a linear function of the main stress. 


6041 ‘between th the United States Air Force, Wright Air Development Division, © 


YIELD CONDITION, 
<4 
7 


and Tt The University of “Michigan, The material sensitive to mean 
stress was completed — in the course of work under | contract DA- -20-018- ORD- 
23276 between the United States Army Ordnance Corps, Land Locomotion Re- 
search Branch, and The University of Michigan, while the writer was a 
_ The contents of this paper formed the basis of a presentation before 
‘Tenth International Congress of Applied Mechanics, in 1960, and an abstract, 


submitted in January, 1960, is to appear in the Proceedings of the Congress. 


APPENDIX, 


‘The following | ‘symbols, for use in the paper, conform ‘essentially 
with “Glossary of Terms and Definitions in Soil Mechanics,” prepared by the 


Committee on Glossary of Terms and Definitions in Soil Mechanics of the Soil 


Mechanics and Foundations Division, Proceedings Paper 1826, October 1958: 
radius. of plate d defined by Eq. 14; 
= cohesive stress computed from. compression test; 
_ cohesive st stress  Compates from extensional test; 


= effective values of cohesive ve stress; a - 
= energy dissipation; 


sy yield function; 


yield moment per unit width | of p plate; 


radius « of plate; 


Lindrical co- ordinate; 


displacemen ment; 


= cylindrical co- 


— 
— 
ar 
— 
= 4 
> 
tm 
on = r/R; 
— 


YIELD CONDITION 


shear stress; 
= angle of friction computed from compression test 


= effective values of angle of friction, 


: 


q 
nd 
— 
| 
— 


— 
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EFFECT OF JOINT ROTATION ON DYNAMICS OF STRUCTURES ab 


By ! Moshe F. Rubinstein! and Walter Cc. Hurty2 


t 


SYNOPSIS 


ss The resuite of a study on the effect of joint rotation on the dynamic charac- 
_ teristics of a multi-story steel-framed structure is analyzed under four levels a 

of assumptions regarding joint rotations. The first ten normal modesand fre- a. 
= quencies resulting from the analysis under each assumption are compared in 


figures and examined 


149 


_ ‘The dynamic analysis of a multi- -story rigid | | frame requires extensive nu- 7 
merical calculations. - Because a multi-story framed structure normally con- 7 
sists of a number of such frames, the computational efforts increase greatly 
when a dynamic analysis of the entire strucutre is required. Therefore, it is 7 
to introduce ‘assumptions: that will r reduce the omout | 


tions, although result in a small sacrifice in accuracy. 
7 ae One such assumption considers the mass of the building — or. concen- 
trated, at the floor levels with three degrees of freedom of motion per | floor 
with translation in two horizontal directions, and a rotation about a vertical 
axis. ‘In the event that the foundation is not rigid enough to fix the i 
against rotation at the base, building may als also rotate about two horizontal 


- Note.—Discussion open — May 1, 1962. To extend the closing date one month, . 
written request must be filed with the Executive Secretary, ASCE. This paper is part 
“of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the 
my Society of Civil Engineers, Vol. 87, No. EM6. 
1 Asst. Prof. of Univ. of California, Los Calif. 
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An additional assumption, introduced by s some investigators,3 considers the - 
girders of the building frames infinitely rigid, thus permitting no joint rotation : 
and restricting the coupling between stories. This assumption simplifies the 
complex spring characteristics of the building frames. However, although lump- 


ing the mass at the floor levels is considered quite reasonable for a a wideclass 
_ of buildings, the assumption of infinitely rigid girders is questioned in the case 4 


_ of multi-story framed structures. The writer’ ’s experience i is that cial 


4 he 


under "such restrictions. The present paper is devoted to a 
A nineteen-story steel-framed building is analyzed under four levels of as- 
sumptions regarding joint rotations. The dynamic analysis is carried c out in 
_ one direction. The building consists of eight frames of three different types 
in that direction. Normal mode shapes and ictal are computed unde under 


1. No. joint rotation takes place (girders are assumed infinitely stiff). 
2. All joints within a floor (for all frames) undergo an equal rotation. 
3. All joints of a given type frame within a floor undergo an equal — — 


al rote 
No restriction on joint rotations. 


‘The first ten modes and frequencies resulting from the analysis under al 


of the above assumptions are ‘compared in — The a analysis was — 
out using the stiffness matrix approach. 


METHOD OF ANALYSIS 


_be into lateral loads, at the floor levels and at the 


“joints. . The v vector of t loads | 5 ( may be related to. the corresponding vector 


of deformations ; through a stiffness matrix 


p 
in es [e] is the stiffness matrix of the frame (see Appendix for the derivation 


of a stiffness matrix for a rigid frame), {yu} is the vector of joint rotations, 


and {y} isthe vector of floor displacements. 


an element fj; of the stiffness matrix [p] has the f following physical signifi- : 
cance. It is the force ‘required at ito permit a unit: deformation at j only ‘fo 


_ other deformation permitted). The terms force and deformation will mean 
: lateral load or moment, and translation or rotation, respectively, as the case 


Vol, 119, 1954, p. 171. 
4 “Forced Vibration and Natural ‘Frequencies of Tall Premes,* od E. 
_ Goldberg, J. L. Bogdanoff, and Z. L. Moh, Bulletin, Seismological Soc. of Amer. , Vol. 
49, January, 1959,p. 33, 


il 
q 
a 
all frame joints with 
— 
— 
a 

al - 
| 

a = 


JOINT ROTATION 


With [p] computed from the geometry « of the frame and the properties of the 
"sections and) 4 ( established from the given loading conditions, the detorma- 


tions may be solved for 


in which [ey is the inverse of the stiffness matrix [p]. 
all moments and shears for the individual frame a 


Using the values of ty ‘ 


members (columns and may be computed. 

> Stiffness Matrix of a Frame with No External Joint Moments. —If the mass 
ot the frame is concentrated at the floor levels and moving to the lateral direc-_ 
tion, only lateral inertial forc forces result, and all inertial moments at the joints 
are zero, or 


‘From Eq. 


— 
7 
— ad 
— 
| 
| 
— 
Substitute this moment value in Eq. 
Eq. 3 can be written in 
{ 0 } = [Rit | { “12 re - 
= [Roi] {9} + [Rea] fy} 
[Ru] = - [Ria] ty! 
| 
‘Substituting Eq. 7 in Eq. 5 yields 
— a} = ( [Reo] - [Rar] [Riz] } 
— 
et 
= 


in which [K] is the modified stiffness matrix of the frame under the assump- 
_ tion of no external moments at the joints. The order of matrix [K] is equal to 7 
_ the number of degrees of freedom in lateral motion, = 
may noted that the of the frame will 1 still “require the 


inversion of ‘the [e] ‘matrix. However, ‘now that [Ri] “has been computed, it 


is worthwhile to take it and use the following relations for the 


| 


= 


: 
Lateral ‘Matrix of a Building.—A relation similarto Eq. 10 may 


“obtained for each frame of a building. Thus, for frames I and II in Fig. 1(b) © 


= 


10, 


ies 
fe 
a 
the symmetrical buildin 
which {x} is the vector of floor deformations in the x direction. Because 
a _ there are two frames of each type, the vector of lateral forces / Hx | for the at 
= by P. B. Morice, Ronald Press, New York, 1959, p. Be 


JOINT ROTATION 
nd Hx Eq. 13 and using Eq. 14 yields 


"7 
x } 


‘The term [K*] i is stiffness for the building in 


The Eigenvalue Problem.—In the building of Fig. the mass concentrated 
at any floorrisdenoted by mp. 
According to D’Alembert’s principle, ‘the: force with a an 


=e my 

Postulating harmonic motion with 


Xp Rp ee 
"Substitute Eq. 20 in 


in matrix 7 


one 


in ‘which a with elements on }on the principal diagonal 
> 


{x} = pm] {x 


premultiply each side of Eq. 23 by [ m m_]-1 


{x} 


a 
EM 6 a 
Substitute }H | 
| 
i ined from the elements of [ Ky | and [Ky | 3 
> 
. 
— f equations becomes, &§ 
— 
g 
: 
— 


= if =e 


(a) Schematic Elevation of ‘Frame 


25 represents the of the classical eigenvalue | problem. For 
a non-trivial solution { x } # {0 io the becenieeauand 4, of the term within the _ 
brackets () in Eq. 27 must vanish 


— {x} = { — 
— 
in which [U] is the unit matrix. 
— 


J For tall multi- story buildings, , Eq. 28 isa high | sii: equation in w, Thus, 
a trial and error solution is resorted to, or the iteration technique6 m may 7 
7 used to solve > Eq. . 2 The non-trivial solution of Eq. 27 yields the natural fre- 
quencies and mode shapes of the building. 
ASSUMPTIONS. REGARDING JOINT 


for the frame ot of Fig. 1(a). Ruy isa square matrix of order 32 equal to the 
- "number « of degrees of freedom in joint rotations. R22 is a square matrix of . 

order 8 equal to the number Of ‘degrees of freedom in horizontal translation. 

R42 and Ro (the coupling matrices) are of order 32- by-8 and 8-by-32, re- 


The to the modified [K] of the frame 


modified stiffness matrix [K] is given all 


no inversion is required. 
Al Joints of a Frame Within a Given Level r Undergo an an ‘Equal Rotation, 
rs Eq. 1 for frame I (Fig. - 


«6 “Dynamics of Framed Structures,” Rogers, John Wiley and Sons, New York, 


1000, 


— q 
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iq. 5 to 
| 
a 
g 
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| fe) 
Beate i | @ 
for 1 = 29, 30,3 


Eq. 1, Eq 
Using Eqs. in in Eq. Ss. 


ji 
1 j= 
Pj, 33 33 * X34 X40 | Pi,40 
iad 
8 
jel 


+ og 
j= 


34 40 
j=5 j-29 | 


Similarly, the la lateral loads Hj, i= 33, 34, . ., 40, become 


= P33 92D, Pag yt -+93 P335+ 
=1 =5 


2 


— if 
— 4 
| 
_ 
— 
3 


4 p* 


: ‘The elements of R11; , Ryo, Roy are identified from Eqs. 34, and are on 
tained by adding blocks of ij coefficients in the original stiffness matrix [ep] 
1). elements of R22 are the same as the elements of 


ord -by- 


_ frames, one matrix of order 8- -by-8 must be inverted in computing the stiff-_ 
‘ness matrix [K] for each frame. A linear combination of the elements kKij of 


“All Joints Within a "Floor (for all frames) an Equal Rotation —An 
equation similar to Eq. 35 can be written for in a building. 


— 


o 
— 
wl 
‘ — 

— 
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‘then Eq. 39 becomes 


[y ]= -4-- is the stiffness is matrix in in the x direction f for the eentire 


building of Fig. 1. The modified ‘stiffness matrix for the entire. under 


( [v2] [van] [v1]? [712] ) 
For the building of Fig. 1, matrices y Y 112,721 and y 99 are each at order 
_* -by-8. Thus only one matrix of order equal tothe number of floors in the 
building» ‘must be inverted in onder to the modified st matrix. 


x Preliminary Considerations .—The nineteen-story building of Figs. 2 and 3 
is a rigid frame steel structure. _ The column and girder sections are indi-— 
‘cated in Fig. 3. The structural properties of the framed members have been | 
- tabulated by the American Institute of Steel Construction (AISC). 7 The floors © 
are constructed of 6 in. concrete slabs. Exterior walls and interior partitions © 
. = lightweight construction. The floor slabs are regarded as infinitely rigid | 
in their own plane so that no girder | change in length takes place. The building 
rests on rigid ground, and the columns are fixed against rotation at the first 
floor level. No foundation settlement occurs. _ The mass of the building is 
at the floor levels, 
The weight ofa typical story Wr = = 1,235, 200 lb. The mass lumped ata typ- a 


Mtypical floor 


—T« Steel Construction,” Amer. Inst. of 
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anine, and the twentieth floor rs 
uivalent to 1.5 times 
Steel Constr., New York, 5th edition, 1947, 
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3 JOINT. ROTATION | 


The dynamic analysis is carried out in the north-south direction. Int this 
_ direction, the building consists of 8 frames (Fig. 2) of three different types. 
Frames on lines 1 and 8 are of one type, frames on lines 3 and ¢ 6 of a second 7 
type, , and frames on lines 2, 4, 5, and 7 of a third type. “—pareaa 


_ The normal modes and frequencies of the building in the north- ‘south di -— 
rection are computed for the followings 
Case a. No joint rotation takes place (girders are assumed infinitely rigid) 
Case b. All joints within a floor (for all frames) undergo an equal rotation. 


Case c. All joints of a given type frame within a floor undergo : an an equal 


rotation, 


Case d. No on joint rotation. 


‘Using Eq. 24 or Eq. 27 ‘the mass matrix fm] of Eq. 43 is is inverted in each © 


7 case. This is, of course, a trivial matter because ['m_] is a diagonal matrix. 


- No additional matrix inversion is required in the computations for Case a = 
matrix of order 19 is inverted in the computations for Case b. ‘Three matrices 
- (one for each frame type) of order 19 are inverted in the computations for — 
 Casec. Two matrices of order 95 and one of order 57 are inverted in the com- 

putations for Case d. This is required because the interior type frames con- 


95 joints each, and the omterice type frames contain 57 joints each. Only 
a In all four cases, the normal modes. and corresponding frequencies : are ob- 
2 tained through the solution of an eigenvalue problem of order 19, equal to the 
_ number of degrees of freedom of motion of the building in the north- south di- 
7 rection. The natural frequencies are are ‘computed by trial and error from Eq. 28. 
For each frequency the corresponding mode shape is computed srom Eq. 27. 
An computations were done with the aid of an IBM 709 computer. _Compu- 


q _ tation time for Case a was 2.39 min; for Case b, 2.64 min; for Case Cc, 2. 69 1 and 


and for Case d, 22.84 min. 


AND CONC LUSIONS 


The first ten mode shapes and frequencies of ‘the building ' computed for 
‘cases a, b, c, and d are compared in Figs. 4 through 13. The mode shapes of 
cases b, , Be: and d are soclose that they are plotted asa single curve in “the: 
figures. ratios are compared in Table 1. 
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JOINT povaTION 
Mode shapes for case a differ from those of case sein. The iad is most 
pronounced in the first, ‘second, and third modes. ‘This is expected, because 
the effect of joint rotation diminishes in the higher modes; it is conceivable 
= the joints in a given floor may not rotate at all in a high mode. Sere tr 
frequencies for case a differ from the corresponding frequencies of 
case d. The deviation is largest in the first natural frequency (109%). — 
= become smaller for the higher frequencies, as the effect of joint rotation | 
‘The frequencies for cases b and c are the same, e, for all practical purp purposes. 


— exterior and interior joints. Once this difference is reconciled by as- 
an rotation for the within a it ‘makes little difference 


This is explained bythe fact that the largest difference in joint behavior is se 


"computed for d are in order of to 6% for the first ten modes, with the 
i deviation decreasing for the higher frequencies. 
TABLE 1. —FREQUENCIES AND ranquenct RATIOS FOR THE STRUCTURE OF _— 


“Frequencies (Radians per 
_Mode | 
Number 


WD 


| 


_ The frequencies computed for ca case a are consistently largest in value; those 
computed for case d are smallest in value; and frequencies of cases b and c lie 
between these two extremes. 7 This is expected, because frequency increases 

7 with stiffness. In case a, the assumption of rigid joints attributes a high stiff- 
“ness to the structure; consequently, high frequencies are obtained. In cases b, 
c, and d, the assumption of rigid joints is progressively relaxed, and, conse- 
quently, the computed frequencies are reduced in value progressively from 
case atobtoc tod(seeTablel). All mode shapes and enaea were com- 
puted to an accuracy of six significant figures, 

In comparing the results of cases a, , b, and c with those of om d, it is con- 
cluded that the assumption of no joint rotation results in mode shapes and fre-_ 
- quencies, with a discrepancy in the order of 100% for the lower frequencies. 
Cases b and c yield very satisfactory results, with only a S%to- 6% discrepancy 
for the ted modes. — Because case b requires less | computations than case c, 


it is 


é 
: 
1 
UCase a se 0 sec "UCase d WwW WwW 
| 3.85 | 1.93. 193 | 1.84 | 2.09 | 105 | 105 
‘| 23.70 12.80 12/10 196 1,06 
45,15 25,99 | 24.86 | 1.82 1.05 
— 51.70 | 31.50 | 31.48 | 30.25 1.71 | 1.04, 
| 
a 
— 
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December, 1961. 
"suggested that the; sonemntann of equal joint ro rotation for all joints swithin a a floor 
be introduced in computing the dynamic characteristics of a multi-story framed 7 
structure. This considerably reduces the computational efforts or computer 


time (a reduction of about 900% in the ‘agen in this paper) with a rather small 


_ computer at t the Western Data ‘Processing Center, University of California, 
‘ Los Angeles, Calif. Without the aid of the computer, ‘this study would not have a 


FOI FORA RIGID ID FRAME 
Consider Fig. 14. Fig. 14(c) is a linear combina ation of Figs. 14(a) and 14(b). 
_ The member ik is moved as a rigid body from position ik to il and then a unit 
=. is applied in the clockwise direction to each end. The stiffness, Kix, 
of member ik at end i is defined as the moment required at end i tc to rotate that 
through a unit rotation with end k restrained against rotation, 
: over factor, Cj,, at end i is defined as the ratio of the moment induced at end 
_ ktothe moment applied at end i , when end i is free torotate and end k is fixed a 
Ka and for end k of member ik are defined in a similar way. 


In deriving the stiffness matrix fora rigid f frame, such as sho shown in 10), 


‘There is no c no change in length of or or girders. 

There is no settlement of frame supports. _ 

- " . Joints are rigid; namely, all ends of members framing into a joint rotate — 
an amount equal to the joint rotation. 


‘The total number of degrees of freedom of motion of the frame is then equal 
_ tothe to total n number of joints and horizontal girder lines (floor levels) in the 


Forces ‘Required to Permit a Unit Rotation at a Frame Joint. —Apply a a 
rotation at joint i = 5 (Fig. 15), with all other joints and floors —_——— 


the definition of ‘Stiffness of a member, over , and an 
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(48), 
Py 5 =- 48) 

- ° Forces. Required to Permit a Unit Lateral Deformation Along a Girder Line.— 

Apply a unit lateral deformation at level r of Fig. 16 with all other floors and 
joints restrained. At end 8 of column5,8 
Mg,5 = - Kg,5 (1+ +Cg5) - 


Similar expressions can ew: written for. in which k = =1, 3,4,6,7,9 (see Fig. 16). 
To obtain in pr + Lr of Fig. 16, equate it to the rth A story shear: rue - 


Mix + Mig 


which the summation extends over columns ik in the rth stor story. 


in which the ‘summation | estente over all columns ik of the (r - 4st t story. an 


To satisfy equilibrium of external forces in Fig. 


lm 
__ To calculate the restraining fc 
Fig. 17. From elementary equilibrium conditions tor columns 
a 4, 
— 
— 
a 
— 
‘ 
— 
= 
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JOINT ROTATION 


aa on Betti’ s Law® to the various load deformation systems used in ibe 


- ii the stiffness matrix, it may be established that 


y computing terms out ofa total of at te terms. 


| 
Th iffness matrix 
= 
— 
4 
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= 
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BEARING OF FLOATING ICE SHEETS? 


‘Closure by y G. G, Meyerhot 


: about the e extent to which the mechanical behavior of ice can be approximated 7 
by a rigid- -plastic material obeying the ‘Tresca yield criterion, Because ice is 
-an elastic-plastic material, it is necessary to distinguish between long-time © 

and short-time stress applications. | For long duration of stressing and slow 7 

- rates of loading, the behavior of ice is similar to that of a ductile material with © 
ee stress in shear of approximately one-half of the yield stress in tension — 
( 


G. G, MEYERHOF, ASCE. —Both H.G, and R. R. “Wood enquire 


12) ,2 thus approximating the Tresca theory of strength = = a 
Under this condition, ice may be represented by a rigid-plastic material — 
_ with an adequate moment-rotation capacity in bending under the limit moment 
Mo, which obtains in practice for long - -time parked vehicles or aircraft and od 
storage loads on floating ice sheets. Although combined loading tests on ice 
_ would be required to check the yield function, the writer believes, from ananal- _ 
_ ogy with other ductile materials, that, at least for long-time loading, the Tresca 
_yield condition applies and is preferable to the simpler Rankine criterion of 
maximum principal stress or moment (“ square” failure condition). Moreover, 
_ previous computations by the writer for concentrated loads at the center and 
_ edge of plates on an elastic foundation have shown that the Rankine failure cri- 
_ terion gives much greater theoretical c collapse loads that are not confirmed by 
observations on ice ‘sheets, The Tresca yield condition was also chosen. with 


a. view to applying the proposed theoretical R saaigele to other related problems 


Under short duration of stressing -and rapid - rates of loading, however, ice 


} “behaves like a brittle material with a shear strength of about one-half of the 


failure, Because the corresponding ‘moment- rotation ‘capacity ‘under 
— yield moment My is limited, the moment redistribution in plates is believed to 
d ibe small and can be neglected in practice for short-time loads on floating ice 
q "sheets. The theoretical collapse load for this condition has, therefore, been 
Fo estimated on the assumption that the ice sheet has either radial cracks or aan 
7 circumferential crack, in the case of concentrated | loads, or a longitudinal © 
-erack, in the case of strip loads. Whereas both assumptions give practically | 
: ‘the same theoretical collapse load for concentrated loads, the actual collapse 


- mechanism favors the assumption of radial cracks, In these cases, the Tresca 7 


October 1960, by G. G. Meyerhof (Proc. Paper 2627). 
ou 4 2 Numerals in parentheses refer to corresponding items in the Appendix Bibliography — 


25 Head, Dept. of Civ. Engrg., Nova Scotia Tech. College, Halifax, N. S., Canada. 
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yield criterion has also been used and would indirect supported by the 
_ good agreement between the proposed — theory and the observations of loading ~ 
aon Gold’s additional information, obtained from field surveys of ice roads and 
landings of pulp and paper companies, is most welcome, and these observations a 
of normal, maximum, and failure loads fall in the same regions as those in-— 
dicatedin Fig. 8. Using the high values of the flexural strength fp deduced from 7 
_ the maximum working loads and the low values of this strength deduced from — 
the failure loads, the various field observations g give a range of fp from about 
100 psi to 300 psi, which is within the customary variation of the strength of 
 fresh- water ice under short-time loads. To the various factors mentioned by - 
Gold in connection with failures of normally adequate ice thicknesses in prac- 
tice may be added the effect of a sudden drop in air temperature and extremely : 
_ cold weather; such weather conditions may lead to cracking of ice sheets due 7 
thermal contraction. . Under such conditions the bearing “capacity of 
‘ice is theoretically about one-half of that of an uncracked sheet (see Fig. 4); a 
_ this would be confirmed by recent field data (22). Before thermal cracks have 
time | to heal, a half load limit should be used on the ice as a precaution. * a 
In view of continued absence of published ‘results of full-scale 
tests on ice to failure, , the series of model tests on floating ice sheets by Gold _ 
7 is particularly valuable. His test results show that the maximum tensile stress 
/ deduced from the yield load at first cracking of the sheets agrees closely with © 
the flexural strength of the ice deduced from the collapse load and the proposed — 
_ equations for central loads. Practically the same flexural strength of the ice 7 
has been deduced from the collapse under edge > and corner loads on sectorial 
4 ‘sheets that failed at a load in almost direct proportion to the center angle of 
the weed and with a crack pattern, as had been suggested in the theoretical - 
a analysis. — . The average observed radius of the circumferential crack of about 
ie: 4Lat failure of the ice sheets is also” consistent with the average value of 
approximately 0.3 L estimated from the proposed theory for narrow wedges. a 
This radius also ‘supports the assumed collapse mechanism of radially - -cracked — 
In this connection, the results of another model loading test are of interest — 
(23) . A floating ice sheet of 8. 4 in, thickness and 26 ft-by-26 ft wide was loaded — 
Bricks the center over an area of 30 in, diameter, the test being made in 5 hr at at 


temperature | of 15° F, The first radial crack» in the sheet was observed at a a 
load of 7,840 lb, when the maximum tensile stress sr is calculated to be 235 
psi (Eq. 1) for a characteristic length L of 200 in, A circumferential crack 
witha radius of about 0,6 L developedat a load of 18,300 Ib, for which the flex- 
ural strength fp of the ice is deduced to be 215 psi (Eq. 38) with a theoretical _ 
crack radius of 0.5 L. Complete collapse of the sheet occurred ata load of 
20, 600 lb, for w which fp is estimated to be 240 psi, The values of fp and Sr are 
vs in good agreement and are of the same order of magnitude as those of Gold’ . 
a ‘tests. The observed deflection pattern is also found to be consistent with the _ 


“theoretical shape (Eq. 45) 
A The influence of the geometry of loading on on the collapse of ice sheets has © 
been referred to only briefly in the paper, and most adhiinns were obtained 
for single circular and strip loads. The analysis for concentrated loads at the © 
edge of semi-infinite plates was, however, derived for a semi- -circular area 
_[ Fig. 5(a)] so that an equivalent radius of V2 a should be used for a circular 


“load of radius a tangent to the edge of an ice sheet before entering Fig. 4(a). 


— 
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DISCUSSION 
"Similarly, wedge plates were analysed fora sector load (Fig. 6) so that an 
equivalent radius of 2 a should be used for a circular load of radius a adjacent 
to the corner (90° wedge) of an ice sheet before entering Fig. 4(a). Moreover, 
_ where the radius of the loaded area is of the same order of magnitude as the 
thickness h of the sheet, an equivalent contact radius at the mid-surface of the 
ice sheet can be used, _ Thus, for thick sheets (a < 2 h) the effective | load | radius — 
_ may be taken as a + h/2, as in the case of punching failure. ae win es 
_ Regarding the bearing capacity of floating ice sheets under combined loads, 
the writer agrees with Gold that under normal loading conditions, the collapse 
load may be estimated from an equivalent circle enclosing the wheel space of 
a single vehicle or aircraft or from an equivalent strip, in the case of a load 
train. This is shown by the comparison between estimated and observed col- 
lapse and working loads of single units in Fig. 8, whereas for the usual contact 
range of a/L between 0. 1 and 0.2, the theoretical radius of the circumferential 


- = FIG. 18.—INFINITE PLATE UNDER COMBINED LOADS 
_ erack of about 0. 5 L to 0.7 L at failure (using the theory proposed for narrow Ww 
_ wedges) compares well with the observed variation of this radius (19). — 7 
ever, where the wheels or tracks of a unit have a wider spacing in relationto _ 
the characteristic length L of the sheet or where a load train is relatively 
short, a more accurate estimate of the collapse load is to be preferred, and the 
- following additional expressions are presented for central loads of various 
common types. The bearing capacity of edge and corner loads may be taken 
. approximately as one- -half and one- quarter, respectively, of that of correspond- — 


Circular Loads. —For two loads of the same radius a and center spacing n 
on an infinite ideally- plastic plate, the combined theoretical flow fieldis shown 7 
in Fig. 18(a ). In the semi- -circular end zones the collapse load may be esti- 
_ mated from Eqs. 4 to 6, whereas in the rectangular central zone of length n, © 
_ the collapse load per unit length may be estimated from Eq. % The corre-— 
sponding total collapse load is, therefore, approximately, 


= (27 +362 


= 
a 
a 
My (fora=0) ......(104) 


_ December, 1961 


= 
sa(1+t) 


or, , approximately, 


fh 
As s the load spacing n oipreninne. Eqs. 104 to 107 approach the upper limit of 
_ the collapse load, which is given by the sum of the collapse loads of the oi 


in Fig. 18(b), and the total collapse load is — 


Mo (for 0.08 <?< ..(109 


Further, the total load of large ice sheets under short- 


(m+n). (1+ +2) 


of the | sum of the loads of the individual circular areas, If 
radius of these areas differs, an average radius a may - used in the above 
_ Rectangular L Loads. .—For a load applied over a short strip of width 2a and 
length 1 n with ‘semi- circular | ends ¢ of radius a. the theoretical flow field at col- | 

_ lapse is identical to that shown in Fig. 18(a) so that the collapse load is, ap- 

- proximately, given by Eqs. 104 to 107. Similarly, for two short strip loads of 
the same size as above and center line spacing m, the theoretical flow field at 
collapse is identical to that shown in Fig. 18(b) so: that the total collapse load 


— =| 3.37 (1+ 3 +L) + Mo ( for 0.05 < = < 1) (105) 
ss Under short-time loading conditions on large floating ice sheets, in practice, ; : a 
a «Eqs. 104 and 105 may be replaced by the following conservative formulas in cs 
accordance with Eqs. 38 to 
ig 
ed flow field is shown 
0) ......(108) 
— 
Age ¢ho land enapings mm and nineroaaca Bac 1028 toa 111 annroaach the wn a 


‘is, approximately, given 108 to 111. As the center spacing in- 
creases, the total collapse load approaches the sum of the collapse loads of the 
individual strips, 
_ Wood, N. C. . Lind, and D. T. Wright — to a influence of deflections on 
‘the collapse loads of plates. Membrane action, change in 1 geometry, and shear- 
“ing deformation of the plates have been neglected in the proposed analysis, be- a 
cause trial by the w riter indicated that their effects would be 
Only where 
are ¢ ne eater than the approximate ‘thickness: of the ice ‘would the theoreticalcol- 
lapse loads be increased significantly, — until supported by field evidence, — 
it seems best to neglect this increase, 
Y. Sebastyan and F. Penner have an extensive | ‘series of loading 
7 tests on concrete pavements of airfields having various subgrades, bases, and 
‘slab thicknesses, Static plate bearing tests have been carried outat free, keyed, 
and dowelled corners and at the center r of pavements when the corners and 
edges were curled down. The test results show that the yieldor first cracking 
loads are much greater than according to the conventional (elastic) method of 
design. For the free corner case, the observed loads agree fairly 


_ the dowelled and keyed corners, the actual set Ai loads are generally greater 
than the theoretical value for edge-loaded slabs. Under central loads, however, 
collapse occurred at considerably greater loads than estimated. ‘This differ- 

ence appears to be largely due to the different yield criterion for concrete and © 
the variable foundation modulus for soils, which may require a modification of : 
the method of analysis and that it is proposed to consider in a subsequent pa- 
per. It is of interest to note that the observed ‘punching loads give a theoretical 
7 tensile strength (Eq. 36) of about one-half of the flexural strength, which is in 
accordance with the usual laboratory test results on concrete. _ Oo 
OA Assur discusses the author’ s assumptions of the position of the yield 
_ hinge circle and the extent of the hydrostatic reaction pressure in relation dl 
some unpublished results of large-scale loading tests on floating ice sheets 
_ These tests show that after collapse the hydrostatic reaction is essentially ll 
fined tothe yield hinge circle in which a circumferential crack occurs near the a 


i location of the estimated maximum ae radial bending moment. As hol 


plastic analysis based < on 1 the actual ‘mechanical properties of the ice. e. In view 
7 of mathematical difficulties, such an analysis is not yet (1961) available. The 

author’s approach is, therefore, based on plastic theory for a rigid plate, which 7 
necessarily implies zero deflection at the yield hinge circle. _ This method, ins 

which an allowance is made for the limited moment- rotation capacity of - 

under short-time loads, furnishes simple solutions for the collapse load — 
- various loading and boundary conditions, and expressions for combined loads 
_ have been added above, as requested by some of the discussers. = = | 

_ Although the actual radius of the yield hinge circle of an ice sheet at col- 

4 lapse is smaller than the radius of the > deflection dish, the latter value has been = 
- adopted in the paper asa safe approximation for the effective radius b in ei 

mating the ultimate load from the proposed theory. If the former radius is 

used, as preferred by Assur, the theoretical collapse loads deduced from hy 
7 3 will be greater, but the difference is relatively small for the usual range of 


contact areas, Further, an assumption that the hydrostatic reaction at ‘collapse 
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extends beyond the yield hinge circle gives | similar estimates, if the effect of 
shearing stresses on the limit moment is taken into account. Assur’ Ss modifi- 
_ cation of Johansson’s theory tends togive excessive estimates for concentrated | 
loads, and for a central point load it gives twice the collapse load estimated | 
7 from the present theory. This difference is due to the use of the Rankine fail- 
ure criterion in the former theory, which does | ‘not apply to ice, as indicated | 
_ previously, On the other hand, the corresponding location of the circumfer- 7 
ential: crack of an ice sheet at collapse is of the same order of magnitude as 
estimated from the equation r= 1.63 va L, based on the author’s s theory of a 
radially-cracked center-loaded sheet. The most useful check of the present 
_ theory may, however, be obtained from a comparison of the estimated collapse © 
loads with the results of loading tests on floating ice sheets. In this respect, — 
the evidence presented in the paper and in the tests by Gold and Linell show © 
_ good agreement between the proposed formulas and observations. A full pub-— 
_ lication of the results of the large-scale loading tests by the laboratory with | 
which Assur is associated would, therefore, be most welcome. 
_ The author is grateful to the various discussers for their kind comments _ 
_ and for adding important experimental evidence; both have helped to assess the 7 
proposed approach to; 2: complex | problem of considerable practical ‘interest. 
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L, ‘Baro HH, Bleich, and P. Weidlinger 


MELVIN BARON, 12 4, A.M. ASCE, HANS H. BLEICH, 13 F. ASCE, ar PAUL 
_WEIDLINGER, 14 F, ASCE.—The last paragraph of A, Eremin’ s discussion 


This statement seems to be based on some misinterpretation of the text by the 
discusser. The paper nowhere makes such a statement, and much smaller in- 
tervals were used. (As explained in the Appendix, stability of the integration 
Pom requires k < T4/n,in which Tq is the shortest period of the beam with lumped 
“masses, . As T4 is much smaller than the fundamental frequency T, an — 


example in lis 


“that may be the. origin of the 

Finally, the writers cannot agree with Eremin’s statement, “. on 
necessary to refer to the period of natural vibration,” because the decay con- 
stant to of the applied pressure must be defined, and its natural (dimensionless) 
definition is the ratio of to the period T. 
‘The writers would like to acknowledge that this paper is the result of inves- 


ee sponsored by The RAND Corporation under contract with Paul Weid- 


a February 1961, by Melvin L. Baron, Hans H. Bleich, and Paul ween site 
Assoc., Paul Weidlinger, Cons. ‘Engr., New York, , and 
of Civil Engrg., Columbia Univ., New York, N. Y. 
3 Dir., Inst. of Flight Structures, Columbia Univ., New York, ‘N.Y. 
Cons. Engr., New York, N. Y. 


DYNAMIC ELASTIC-PLASTIC ANALYSIS OF STRUCTURES? | 
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ee ; s mericalexample it is stated that the decay constant of the pressure wave inthe ‘i 
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CESAR A. SCIAMMARELLAI6 and AUGUST DURELLI.17—Mors 


on areas of comment: 
* "Objective of the paper (first paragraph); 


Limitations of the method (second paragraph) ; 


Each point will be considered in turn. 
Apparently, Morse has missed the objective of the paper. Noting that 
the mathematical theory of deformations has its starting point in the displace- a 
ment field, it was the objective | of the paper to present Moire data in terms of 
this ‘displacement field and so allow the computation of the deformations from — 
moire data. No “new method of measurement ... with greater  seeeealll was | 


8 measurements he describes ‘fifth, , and sixth 


2 Both large displacement and small displacement equations were pre- 
sented in the paper, and it is the judgment of the writers that there are many 
cases in which the small displacement equations may be used, The question of 
what amount of strain is necessary to produce fringes depends on the line den- 4 
sity of the available grid, the geometry tobe analyzed, and the maximum allow- 
able strain for the material used. Further examples of problems with large 
strains and large rotations are needed before definite conclusions may be drawn) 
as to whether Eq. 15 or 16 should be 1 used, rather than Eqs. 20 and 21. - iit 
3 _ The distinction between the Eulerian and Lagrangian descriptions ismade _ 
in the section “Fundamental Properties | of the Moire Fringes” and the ‘meee 
following. 
a 4, The correction of Eq. 19, although not engeg any of the ‘conclusions, _ 


help in the understanding of the presentation 


i In comparing the “geometric | approach” y with t the use of Eqs. 15 and 16, a 
it should first be pointed out that the field of u and v , values completely — 4 
the strain field, and that the fieldof 8u/ax, du/dy, av/ax, av/ay also defines 
the strain field; thus, the four other measurements Morse describes are — 


a cima 1961, by C. A. Sciammarella and A. J. Durelli (Proc. Paper 2736). 
16 Argentine Atomic Energy Comm., Buenes Aires, Argentina, 


7 Prof., Civ. Engrg., Catholic Univ. of America, Washington, D 


MOIRE FRINGES AS A MEANS OF ANALYZING STRAINS@ 
Closure by C. A. Sciammarella and A, J. Durelli 
— e has 
Use OF am, Lapranpian tsecond paragraph) ; 
| 
q 
x 
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sufficient to construct the triangle. wee 
_ The one advantage of the geometric method, cited by Morse, : is the use of - 


one photograph to determine some | values of ‘strain, Th This advantage results 
al ‘all a uniform strain field at the point of analysis. In order to analyze strains 
with one photo, it must be further assumed that rotations are negligible. © —Ob- 
- viously, if there is rotation,a shearing strain along the model grid lines would 
not affect the fringe pattern, but would change the value of strain perpendicular 
t the master grid. In other words, the u-field does not t uniquely | determine — 
if there is ‘rotation, Cte 
us, it is is obvious th that regardless a the e approach, g generally, we 
J —— fields are necessary to compute strains, . In the special case of no rota- an 
“tions, as on . lines” of symmetry, either approach m may be used to determine | 
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str rains from one photograp 
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WIND. INDUCED VIBRATIONS IN ANTE A MEMBERS@ 


‘Closure by William Weaver, Jr. 


WILLIAM WEAVER, ASCE. — The writer agr agrees with Woodruff that 
there appears tobe no “ critical” velocity above a ney ae number of approxi-_ 
mately 3 by 105, 

read The values of Cx reported inthe experiments we: were obtained by direct force 
measurements for stationary cylinders. Values of Ck for oscillating cylinders — 

a and flexural members were calculated using the incremental load factors(ILF) _ 


‘defined in detail in the paper. The | value. of 0.20 for large stacks anes by 


‘Woodruff corresponds fairly well with the writer’ sresults, | 


With regard to Woodruff’s remarks concerning the frequency of oscillation, 


all responses of cylinders and flexural members were observed tobe at natural 


- question was raised ont the numerical value of ‘R in Appendix i, but the 


which, R refers to the number; V is ‘the velocity, in n per sec- 


ond; D denotes the diameter, in feet; vis the velocity, in miles per hour; a 
describes the diameter, in inches; and v = is the kinematic viscosity of air, in 


square feet per 


Pietra 1961, by William Weaver, Jr. (Proc. Paper 2748). 
4 Asst. , Pet. of Civ, Engrg., Stanford Univ., ‘Stanford, Calif. 
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RESPONSE OF MULTI STORY ORY STRUCTUF ES 


Discussion by John A, Blume 


“JOHN A A. Bi BLUME, F, AS F, ASCE. — Work in the detailed of inelastic 
-storied buildings to earthquake ‘motion should 
Berg, F. ASCE, is to be commended for his efforts in this ett field. The 
- problem is one of many parameters, but carefully controlled research may © 


lead to the isolation and evaluation of the basic factors, Although many of these 
a parameters have been considered by other investigators, and analysis tech- 
niques have been proposed to include the effects of damping, multi- -story be- 
_ havior, inelastic response, and the combined participation of frames, shear 
walls, and other building elements, ‘more data are needed to refute, substanti- - 
ate, or monny the proposals and to provide a needed background of reference 
‘The fact that most structures respond to major earthquakes inthe inelastic 
range and that this response should be considered either directly or implicitly 
in design practice is gradually being more generally recognized, | There still - 
is a great tendency, however, for many designers to assume that damping, or 
perhaps a few shear walls, will take care of all the problems. The numerical 
evaluation of these and other important factors is thus desirable not only for 


basic but also to help overcome inertia in the practice of building: 


consider the implications of the simplifications and to compare the model char - 


acteristics with those of actual buildings. The zero damping situation does 

not, of course, occur in structures, However, it has been found5,6,11,16 that 
damping at small amplitudes is not great, and values of 5% to 10% of critical 
damping have been proposed as applicable to the non-damaged state. 6,17 Al- 
though greater energy dissipation that that represented by the 5% to 10% vis 

- os cous damping occurs in the inelastic range, it seems desirable to differentiate 


‘The author’s are idealized, as they should be. It may be 


without damage to or ‘permanent ‘distortion of the structure, and ‘the inelastic 
range of cracking, yielding, or other damage. Berg has used viscous ee 


_ conditions that result in n 14. 8% of critical for the first mode of the four-story © 


4 April 1961, by Glenn V. Berg (Proc. Paper 2790). ethos. * 
18 Pres. ., John A. Blume & Assocs., Engrs., San Francisco, Calif. Wea 
_ 16 «4 Structur al-Dynamic Investigation of Fifteen School Buildings Subjected to Sim-_ 
ulated Earthquake Motion,” by John A. Blume, ‘Sharpe, and Elsesser, prepared for the 
State of California, Dept. of Public Works, Div. of Architecture, distributed by Printing: 
Div., Documents Sect., Sacramento 14, Calif. 
17 “A Reserve Energy Technique for the Design and Rating of Structures in the - 
elastic Range,” by John A, Blume, Proceedings, | Conf. on] Ear thquake Engrg., 


Tokyo, Japan, 1960, 
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7 _ model and the other 6.25% of critical for the eight-story model. Because _ 
nominal value of 5% seems applicable to the frame and from 5% to 10% for the 

‘shear walls, depending on detailed conditions, the values used are not —— 

_ cantly greater than those often considered applicable. However, it would seem — 

well to include a slightly damped frame among the elasto-plastic models. | 

- The use of zero damping values is necessary to compare the author’s re- 
"sponse of the elastic and inelastic frames. This is unfortunate, not only be- 

- cause zero damping is not directly applicable to structures, but also because 
undamped response to earthquake motion may lead to much greater peaks and 
valleys and irregularities than for a more typical value such as 5% of. critical, 7 

that “irons out” many of the irregularities of response. 
_ Berg notes fundamental periods of 2.01 sec for the four- story frame and _ : 
1.45 sec for the same frame with the added shear walls. These are extremely 
long periods for any four-story building and, in fact, are longer than the tradi-_ 
tional type fifteen-story building, studied by the writer, and longer than a mod-— oe 

ern fifteen-story building, for which periods were carefully determined all 
- _ through construction, 18 © The height squared divided by the width of the four- } 

‘story building is approximately 45, Building periods in this category19 of ) 
_height-to-width generally range from 0. 3 sec 0.8 8 sec and and usually would 

Anattempt was made by the writer to periods from the data 
shown in Figs. 1 and 2. Assuming all columns oriented with their greatest mo- 

‘ments of inertia in thedirection parallel to the bents illustrated, and neglecting ~ 
all resistance but the steel frame, the following fundamental periods were ob- 
tained, using a a portal method and the procedure: 


ous! —_ 3 our-story frame wi ips per in, shear shear = 


wall stiffness added to each each story 1.24 sec 


_Eight- -story frame 84 sec 


The eight- story periodalso seems very long in comparison with average build- 
ings. It would be well for Berg to list his shear spring constants and elastic © 
limit shear values so that readers would have available the same dynamic ic pro- 


perties that were used inthe experiments. 


sponding stories with fireproofing, stairways, and other. ‘elements 
_ cluded, this does not detract from the value of the response data obtained for 7 
multi-mass dynamic systems of | specific characteristics . It does indicate, 


7 however, that the results should not be considered typical for four and eight- 4 _ 
story buildings. The models might be “extended” to approximate dynamically " 4 
equivalent buildings of of more stories. “procedure of coupling stories was 


- “Periods o ofa a ‘Modern Multi- -Story Office Building During Construction,” by John ¥ 
A. Blume and R. W. Binder, Proceedings, Second World Conf. on Earthquake Engrg. a i 
Japan, July, 1960.00 | 

19 “Lateral Forces of Earthquake and Wind,” by Arthur W. Anderson, John A. Blume, 
Henry J. Degenkolb, Harold B. Hammill, Edward M. Knapik, Henry L. Marchand, Henry 
_C. Powers, John E. Rinne, George A. Sedgwick, and Harold O. Sjoberg, Transactions, 

“ Earthquake Stresses in Shear Buildings, discussion by John A. ‘Blume, Trans- 


actions, ASCE, Vol. 119, 1954, 194, 
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DISCUS 
used ina a ‘research effort6 to reduce a | building toa 
mass equivalent system that could be accommodated in the available electric 
analog. The writer has taken the liberty of doing the opposite with the author’s — 
i -story ‘frame, as shown in . Fig. ii, where, as indicated a’ ‘the dash lines, it 7 


EXTENDED 


Fig. 11 shows elastic acceleration cosfficients teow two sources for the 1940 


Centro (N- s) earthquake, Only “all-mode” results are Plotted. The natural 


the shapes a and values of the various elastic ae, when all factors are 7 


be considered in making comparisons, ~ “There is good reconciliation between 


| 
— 
Jy | | | | 
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analyze and compare all available purty A few comparisons will be made herein: 
a ya The upper three stories of the four and eight-story framed structures are 
> 


identic al in weight and dimensions and vary only slightly in column nn properties. — 
. The response of these stories, as shown in Table 1 for the Helena a earthquake, — 
is considerably greater for the eight- story undamped model than for the four- | 
model; this is as would be expected. It is noted, the 


(1) ELASTIC. FRAME. 
ELASTO- PLASTIC FRAME 


ELASTO- PLASTIC FRAME 
o TOTAL DEFL.OF | 


PLASTIC FRAME 
om 
ELAPSED TIME FROM 


START OF EARTHQUAKE 

ABOUT 4 ‘SECONDS 


~ 


FIG. 12. _—FOUR -STORY UNDAMPED FRAMES (RESPONSE TO 


EL CENTRO MOTION FIRST YIELDING) 


"frame of the four-st -story model (damped to 14.8% of critical, whereas the eight- 


story is only damped to 6.25%) has essentially equal or greater seismic coef- 


ficients, in spite of lesser deflections . These particular damped results | seem 


The elasto-plastic response in Berg’ s “paper is limited to the four- -story | 
models under the 1940 El Centro N-S ground motion. Col. 4 of Table 1 gives 


21 Discussion by John A. Blume of “Behavior of Structures During Earthquakes, > 
ASCE, Vol. 86, No. EM 3, June, 1960. 


22 “Dynamic Response of Elasto-Plastic Frames,” by Joseph Penzien, 
ASCE, ‘Vol. 86, Ho. ST July, 1960, 


tm 
considered, It isto be noted, however, that me elastic accelerations for models 
of shorter period would be considerably greater,16,17,20,21,22 
It isimportant to obtain response data for all possible multi-mass systems, 
: 
13 
| 
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DISCUSSION 
the maximum deflections and accelerations for the undamped e lasto- -plastic 7 
frames; these values are less than for the undamped completely elastic frame, 
as shown in the second column. _ This is interesting when questioning13,17 


elasto-plastic response phenomena. In this case, as the author notes, the max- 
imum deflections under elasto-plastic response are greater than for the cor- 
earthquake exposure of the elastic system. However, 


“the elasto-plastic yieldings “may be confused by cumulative effects. In order 
to explore this problem and alsoto attempt to evaluate the relative energy de-_ 
mands, the writer has scaled values from Figs. 9 and 10 for one instant of time 
and for all stories. The p peak va values occurring at about 4 sec elapsed time — 
were selected, because this apparently represents the first excursion beyond | 
the yield point and is thus uncomplicated by prior lurches and any cumulative 
— Although no great accuracy was possible o: or is claimed for Fig. 12 and the 
related results, it is believed that the general trend is apparent. It is noted a 
- that the completely elastic ‘system shown as curve (1) consistently has less 
- total deflection than the elasto-plastic system of curve (4). Curve (3), that is 
the deflection beyond the yield point for the elasto- -plastic frame, is interesting 


‘fourth | stories had m minor ‘inelastic 1 response 2 as compared to the first and thir 
Fig. 13 indicates the shear and deformation relationships obtained for the 
same (4-sec elapsed time) instant. Obviously, the work done in the inelastic 
: range, for this case at least, varies considerably between stories. The term — 
C' represents the seismic coefficient at which yield occurs as obtained by this 
Scaling procedure, with the aid of the elastic spring constants represented by 
the slopes in the figure. _ The low yield value of the third story seems incon- 
‘sistent with the other stories and my be due to the difficulty of scaling or the. 
model characteristics, 
‘The dash line in Fig. 13 roughly snisaiain a mode or distortion shape, be- = 
cause the deflection scale is the same for all stories. This appears to be a 
fourth mode of vibration, whereas a line connecting the yield points would in-| 
_ dicate a third mode, Such modes are , of course, reinforcing the fundamental 
‘mode motion. _ These conjectures are not too unreasonable, considering the 
facts that the spectral peak of the El Centro earthquake is in the period range 
- of 0.2 sec te to 0.3 sec, the third or fourth mode of this flexible model would be 
in the same general spectural region, and Figs. 9 and 10 plainly indicate short 
period motion, as well as the 2 sec fundamental-mode motion, = = ———— 
_ Story strain energy and inelastic work were computed by the writer for the 
: i conditions shown in Figs. 12 and 13; that is, for ‘the instantaneous deflected 
- position. In addition, strain energy for the elastic response of Fig. 9 was com- 
puted for the same instant; ; results are shown in Table 2. <1) tag 
‘Table’ 2 indicates that there is somewhat more energy in the elasto- plastic. 
"case when it is not completely elastic, especially in the lower story. The dif- 
ferences are not great, however, in view of the scaling and perhaps other in- | 
7 Moreover, energy ‘is assumed to affect design a as a square foot 


rather than as a direct function. ‘Sa 
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4 The reserve -energy technique6,17,23 for inelastic analysis and design re- 
duces elastic ‘spectral accelerations to. equivalent inelastic values in 
code terms, with allowance for energy absorption during excursions into the 
inelastic range. For an ideal elasto-plastic case, , such as Berg has assumed, ~ 
value of R, the reduction factor, for a single mass system becom becomes simply 


FIG. 13. —STORY SHEAR AND DEFORMATION AT FIRST | 


1/2 u-1, which is the total deformation divided by the yield point 
hy oo Such values are shown in Table 3 for the data obtained for the same | 


4 

23 “Electrical Analog for Yield Spectra,’ wi discussion by . John A. Blume, 
Proceedings, ASCE, Vol. 86, No. EM 3, June, 1960. ah 


A 4THSTORY | 
= 
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"DISCUSSION 
et The tabulations in Cols. 6 and 7 of Table 3 are made u under the wesiiiie 

_that each story is adiscrete or independent dynamic system, that is, of course, -- 

_ not the case, except in the fourth mode, in which this approximation may have 
some merit. It is interesting to note, however, that the values of Col. 6 are 
not much different than those of Col. 4 and, also, that the Col. 6 energy values’ 
fall closer to the test data of Col. 4 at all stories than do the values in Col. TS 

that ar are based on deformation, | rather than energy preservation, in the inelastic | 
range. This concept of stories being treated as independent inelastic systems 

deserves further exploration, It has been ‘suggested!? for first stories, ome 
aoe conditions, but not for other stories, It would ‘seem that, generally, for 


Berg’s results, ,energy preservation more closely represents the test data than — 


dis lacement reservation, 


“TABLE 2 ew DEVELOPMENT, IN N KIP-INC HES, AT 4 SEC ELAP SED TIME 


Elastic Frame 


| 8 
115 


TABLE 3. ACCELERATIONS AND REDUCTION FACTORS 


Experimental Data |. Calculated Values 
Elastic ‘Yield 
Acceleration ‘Value, 
— 
a, in 
gunits cr 


a 


under ‘consideration. Is the total energy ‘capacity above each story effective, 

_ or should it be reduced because of the improbability of all above stories par- 
_ticipating at full capacity ? Is the value of stories below the one under consid-— 
- eration also effective and, if so, in what t degree ? ? There is no simple answer, 

7 secre since time- -history of performance isa factor, unless, perhaps, the — 
7 = story consideration noted previously might be found to have practical 


application under conditions. The reserve - energy technique 


if 
ey 
45 
o 
67 
165 
— | 
j 
1 | oz | on | ov | is | o6 | ose 
basic problem in multi-story building considerations inthe inelastic range ©: 
oe is how much the energy dissipation in other stories helps or affects each story a 
ig 
= 
 &g accordance with various probability factors. The assumption of a straight line i 
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deflection but probability reduction for | energy, produced a first-_ 
story R factor of 0.78 in this case, as compared to the 0.75 found. The R- 


factors for the ‘upper stories were found to greater than the ‘experimental, 


is mow that extensive research is) needed to provide complete ai answers 


_ be correlated with actual field conditions and structures, considerably extended, 
and the data ata should be analyzed in great detail. The writer’ S approximate anal- — 
_ yses ; presented herein are not intended to offer or to imply any sweeping con-_ 7 
clusions from this limited work but to suggest types of detailed considerations, © 7 
_ taken step by step, that seem to be needed in this problem. _ Unfortunately, in 7 : 
- the current state of progress, with different investigators using so many dif- a 
ferent parameters, comparison and reconciliation cannot be made effectively 
For example, Berg used four and eight story models of specific periods and — 
_ the north-south component of the 1940 El Centro record, whereas Penzien22 
used the east-west component and six-story models, A coordinated program. 
‘ could, no doubt, lead to earlier isolation and evaluation of the many parameters, 7 


done the author, by Joseph Penzien, 22 in New Zealand,24 and in Japan 


recognized that there are several types of damping ina building and 
that truly 3 viscous damping may account for but a small part of the totalenergy 
absorption. However, it seems desirable to continue with the mathematically © 
_ convenient concept of viscous damping, but to confine the values to. the low | 
~ amounts found at small amplitudes (perhaps 5% to 10% of critical) and to allow a 
_ for greater energy consumption by friction, grinding, and so forth, in the dam- 
aging or yielding range, as another parameter, that is, work done, This is 


has been done. It is not ‘necessary to ust use equivalent viscous 
Bein implies, to analyze structures with the aid of response spectra. The 
reserve energy technique 5,16,17,23 may be used with any spectrum desired 
and any viscous damping value, although 5% of critical is generally recom- 
mended, and convenient smoothed spectra are provided. _ Although this tech- 
nique needs refinement and more test data before its generaluse or its use as 
= exclusive design procedure on multi-story buildings would be considered, 
itis general in application, yet quite simple and adaptable, and has been used a 
_ extensively in the analysis of all types of buildings. One of its benefits is that - 
it separates damping, damage or yield, and | deterioration. ‘By combining these 
into some such broad term as ; equivalent damping, identification is lost, and 
designers lose the art or “the feeling” of design, 
‘There may be little residual doubt of, but perhaps some reluctance to ac- 
cept, the fact that most multi- -story buildings will have to enter the inelastic 
_ ‘range to resist major earthquakes. The principal problem is, and this has © 
been developed extensively,6 that modern buildings without the aid of “non- — 
walls, partitions, concrete fireproofing, and so forth ‘cannot be ex- 
pected to have the reserve energy capacity of the traditional type buildings or 


more general acceptance of the f facts. of structural dynamics and more realistic 


_ design analyses, at least for special cases. Deformation and energy are basic 


8 perform so well as they did. Not only is more research needed, but also” 


24 “Electrical , Analog for Earthquake Yield Spectra,” by G. N. Bycroft, M. J. Murph ; 


as J. Brown, Proceedings, ASCE, Vol. 85, No. EM 5, October, 1959. 


y 
— 
fa 
. 
— 
| 
Bee 
— 
q 
= 
lm 
‘ 
— 
il 


rms 
not only of primary loads but also of secondary effects, eccentricities, erm 
tendencies under severe story deformations. The interaction the 


but when contemporary architectural designdeletes most of these elements and 
materials from the buildings, the design professions are faced with a serious 
problem. Modern buildings require more advanced earthquake designs than 
their predecessors, if they are to survive severe earthquakes, The codes and 
code coefficients along may not be expected, nor are they intended, to —- 
‘sede good engineering judgment based upon all available and pertinent facts. 7 


This somewhat extensive discussion is not intended to be critical of this © 
excellent research effort but to emphasize certain factors considered impor- 
tant in the basic, but little- recognized, subject of inelastic response of multi-— 
"story buildings to earthquake motion, to ) encourage more efforts in this field, 


terested n may ¥ work for the or the public welfare. are. 
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‘CONTINUOUS BEAM-COLUMNS ON ELASTIC 


SENG LIP LEE,12 M, ASCE, T. M. WANG,13 and J. s. KAO, 14_ The three- | 
moment equation for continuous beam-columns on elastic foundation derived : . 

- by Lu, Eq. 60, defines the relationship between the end moments and deflections. 

In comparison, the use of the slope-deflection equations, Eqs. 18a and 18b or 

Eq. 22, as the case may be, yields the relationship between the end slopes and 
deflections. The former leads, in general, to a system of simultaneous equa- | 


tions that involves more unknowns than the latter. For the problem treated in 
Example 1, for instance, , Eqs. 62, 63, and 64 contain three unknowns Ma, Mb, — 
and Yb» , whereas Eqs. 30a and 30b involve only two unknowns @p and Yb. ye 
‘The proposed method | may be ‘extended to the analysis ‘of continuous plane 
frames with members which are elastically restrained. In this case, the use - = 
the slope-deflection equations generally introduces three unknown displace~_ 
ment components per joint, the rotation and two deflection components. On the 
other hand, the application of the four-moment equation, which may be derived : 
similar fashion as Eq. 60, involves more unknowns per with ex- 
— ception of joints with only two members. 
five unknowns, three moments, and two deflection the 
_ fourth moment being expressed a priori in terms of the other three. On the 7 
| Sandpoint of the number of unknowns involved, therefore, it must be concluded - f 
that the use of the’ slope - -deflection equations is preferable. ews 
_ As stated at the end of the conclusions of the paper, the moment-distribution 
-method should be used when large number of joints are involved. In this re- 
spect, the stiffness influence coefficients contained inthe slope- deflection equa- 
tions furnish the distribution and carry-over factors, whereas the fixed-end 
_ moments give the residuals in the joints to be balanced by successive approx 
discussion takes nothing away from the significance of Eq. 60, 
- which gives the general form of the three-moment equation and should be use- _ 
_ ful in cases in which it is desirable, for some practical reasons, to formulate 
‘the problem in terms of the end moments, instead of the 
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a April 1961, by Seng-Lip Lee, T. M. Wang, and J. S. Kao (Proc. Paper 2801). 

12 Prof. of Civ. Engrg., Northwestern Univ., Evanston, Il. =~ oe 

13 Asst. Prof. of Civ. Engrg., Univ. of New Hampshire, Durham, N.H. 

Walter P. Murphy Research Fellow, Civ. Engrg. Northwestern , Evans- 
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_ JAMES A, LIGGETT, 8 M. ASCE. - —The writer finds Wasley’ s division of 
zones I and II rather strange. The disturbance created by the end of the c curb_ 
_ does not only travel in the transverse direction but in all directions. It is -— 
this reason that the author was forced to restrict his solution to supercritical — 
flows. Therefore, the equation introducing the | small wave velocity Eq. 6 ‘should 


= Voz* 


the author is trying to find the “equation of the characteristic 
(line along which small disturbances propagate) that begins at the end of the 


. The equation should be 


“for v= 0. a Eq. 20 follows directly tvem the definition of the small wave angle? 


which F is the lc e local Froude number in this F 


; of the author’ 's , development will show that he has neglected the c in the rian 
-inator of Eq. 20. Thus, the end of the curb will affect besa flow upstream of ‘the 


line that he has indicated. BS 
—— Eqs. 2 and 19 to integrate Eq. 20, one ‘abies 
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Paper 2880), 
Asst. Prof. of Civ. Engrg. Covell Univ., Ithaca, N. Y. 
ge Elementary Fluid Mechanics, ” by J. K. Vennard, John w jiley Sons, Inc., ‘New 


York, 4th p. 


HYDRODYNAMICS OF FLOW INTO CURB-OPENING INLETS@ 
“Discussion by James A. Liggett 
— 
q 
19) 
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A short analysis 
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which would rey replace author’ s Eq. 10. "Accepting the author’s 
Eqs. 13 through 1| 16, Eas. 17 and 18 would — 


| 
y 
+ 
“The corresponding equations in the second column of > Table 1 (under “Dimen- 
Bees ) would also change, but the equations in the third column remain the —_ 


same. Hence, the author’s non-dimensional plot < Fig. 9) is correct 1 for the 
‘theoretical points, b but the experimental ‘points would be plotted at _ a 4 


_ instead of using Eq. 18a as the normalizing function in the x- direction, 7 7 
oe: It would be interesting tosee if the present analysis brings the experimental 


work into closer agreement with | the theory. The change seems to move the 


plotted experimental points outward on the R- axis. 
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by Campbell Masse 


CAMPBELL MASSEY. 5— analogue suggested by Freiberg for obtaining 
influence lines of moment, shear, rotation, and deflection in beams is much 
strips, are ‘variable rather than fixed in value. A beam varying in flexural 
‘rigidity along its length may then be treated by the analogue simply by setting 
‘the values of R' proportional to (E I) /(ds) of the relevant structural element. 
“ail It may | be of interest to report that an analogue similar to that described in 

the ‘paper was designed and built in the Civil Engineering Department of the 
‘University of Western Australia during 1958 and described in a paper presented 
at th the Congress of the Australianand New Zealand Association for the Advance- 
ment of Science, held in Adelaide, South Australia, in August of that year. Full 
_ details were subsequently presented by R. Leicester, 6 whose paper was tad 
viewed in Applied Mechanics” Reviews? and wh who extended the ie argument to in-— 


_ The analogue is used for undergarduate instruction in the University of 
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a August 1961, by S. Freiberg (Proc. Paper per 2909), —— 
5 Senior Lecturer in Civ. Engrg., Univ. of W astern, — alia. _ 
6 “Resistance Analogue of Beams and Arches,” by R. Leicester, Journal, Inst. Engrs. 
of Vol, 32, No. 4-5, April-May, 1960. 
Review 1401, Applied Mechanics Reviews, | Vol. 14, No, March, 1961. 
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